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ON THE REAL ROOTS OF A RANDOM ALGEBRAIC POLYNOMIAL 


S. BAGH 
Department of Statistics, Sambalpur University, Jyoti Vihar, Burla 768019 


(Received 3 May 1988; after revision 6 September 1988) 


et ¥.) 


neg De asequence of identically distributed independent random 


variables belonging to the domain of attraction of the symmetric stable 
law (Gnedenko and Kolmogorov, p. 171). Suppose ao, a1, --- @n are non- 


zero real numbers and max |lar{=k,. min |ar| = t, and aes O 
0Srsn O0<rSn tn 


(log n). If Nn (@) be the number of real roots of the equation z UeAGX Ue 
r=0 
then for all n > ng, it is proved that 


log n 


32 log (( ) dog n)al® ) } 
n 
161 e log (( es ) (log nsi* ) 
S : 


log n 





P ie Nn () < 


< 


INTRODUCTION 


§1. The object of this paper is to find the lower bound of Nn (w) of the random 
polynomial 


t (x) sap ay Xy (~) x" err} 


when the coefficients are not identically distributed and belong to the domain of attrac- 
tion of the symmetric stable law with index ‘a’ i.e.0 < «<2. The problem of 
finding bounds for Nn () has been considered by various authors. Such type of equa- 
tion has been considered by Dunnage*®. Samal and Mishra® have considered the same 
case when the variance is infinite. Mishra, et a/.6 considered the lower bound for 
Nn() when the coefficients X;’s belong to the domain of attraction of the normal 
law. In this paper we will extend the results of Mishra ef al.6 to the case of domain 
of attraction of the symmetric stable law. 


n 
Theorem—Let f (x) = % ar Xr x" be a polynomial of degree n where Xr’s are 
r=0 


identically distributed independent random variables belonging to the domain of 
attraction of the symmetric stable law with characteristic function exp {— | t | * 4 (#)}, 
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0 < « <2,h(t) being a positive slowly varying function in the neighbourhood of the 
origin. Let ao, a1, 42,... @n be non-zero real numbers. Then there exists a positive 
integer 0 such that for alln > no, 


0 
p 5 : Nn (w) < aie (( i) (log =) 


k 
161 e log (( ie ) dog ny) 


log n 
kn ) rao = , 
im — ]isfinite,wherekn = max |a,|andtn = min |ar|. 
poe ae ( In O<rSn O<rSn 


§2. For the proof of the theorem we need the following definitions, notations 
and lemmas in the sequel. 


Let M be the integer defined by 


M = ( = )s (log n)s | +1 eels) 
n 
where [x] denotes the greatest integer < x; and let k be determined by 
M2 cn & M242, SAG) 


From (2) and (3) it follows that 





bls 5 Eid 5) 13 as ee ee —__¥8. (eee 
4 log (( =) (log n)5/* ) 2 log ((—* ) (log nyo 
n 
We consider f(x) = 2 ay Xy x" at the points 
r=0 
1/a 
tha decal) FAP M2m at) 


k k 
for m =f Jel Ss jaa. 


There are k/2 points if k is even and (k + 1)/2 points if k is odd. We express f (xm) 
for sufficiently large n as sum of the three parts as follows: 


J (xm) =(Z + F+ B )ay X (w) x" 
1 2 3 id 


where the index r ranges from M2m-1 4 | to M2m+1 ind 


7 from 0 to M2™-1lin D and 
2 


from M2m*1 + | to nin 3S, 
3 
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We write 
Am (o) = z ay Xp (w) x” ...(6) 


and 


Rm (@) = ee + 2) dy Xp (w) x", ute) 


Obviously Am and Am,1 are independent random variables. 


Lemma 2.1—If h (t) is a positive slowly varying function in the neighbourhood 
of the origin, then for p > 0, 


(i) lim te h(t) = 0 
0 


t> 
(ii) lim t7*h(t) = oc, 
t>0 
These results follow from Karamata’s theorem (cf. Ibragimov and Linnik?, p. 795). 


We define normalizing constants Vm by 


V~ = |ar|* se h (ay a 6/V'm) 
1 


m 


where @ is an arbitrary small positive number. 
Lemma 2.2—Vm > tn (2 mm J ford > 0. 


Proor : A(t) may be bounded or unbounded whens — 0. If lim A(t) = ©9, 


t>0 
then there exists ; > O such thathA(t) > 1, fort< 1%. If (t) is bounded, as A (t) 
is positive in the neighbourhood of the origin, there exists d > O such that A(t) > d. 
Hence A (t) > d for both the cases. So for large n, we get 


oe ——__,—__., ...(8) 
log n log (( ae ) dog nj) 


Now 


gg Rg “ lar {> SP (ay x 6/Vm) 


m 


2m41 
ar .& 
> d Xm ‘. 
M214 
mz” 
D } ar 
> d xm 6 
mem-1 


(equation continued on p. 658) 
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| 
a3 dt® (M2" — ment (1 a= qx en 


d « 9 
— am 
> - a Me, 


Thus Vm > tn ( : M2m jem, 


#9) 
Lemma 2.3—Let 
2 V 
h= fo: | Bar X (0) x, fees et 
3 
V: 
T2 = jo: | Bae X; () ee [ oe 7 } 
T= {o: | Rm() | > Vm} 
and G be the set of all points for m = | a + 1, Ee + 2, ... k where 
| Rm (w) | > Vn. 


Then 


ot (( A) dorm) 
Pi(G) 129 6 
log n 


Proor : The characteristic function of = a, X;(w) x’, is given by 
3 


pm (t) = exp {— | t | * hm (1)} 


where 
ef « ar 
Pa (C) Ree bd ( arx},1 ). 
Now 
4[Vmn 
Vn 
P(™M™) <2— | Ye | dm (t) dt| 
—41V 
(cf. Gnedenko and Kolmogorov, p. 54) 
4IV in 
V. 
< fli - sme iat 
—4/V¥ 


But 
hm (t) = - | ar | * x" h (ay x" 1) 


(equation continued on p. 659) 
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<|tl-eZ | a, | ae x! 0) 
. m Ne 


Since by Lemma 2.1, for P > 0, 


Fr 
h(arx,, th< | ar x t|-* ast +0. 


But 
(1 ~ gm(t)) = 1 ~ exp(— | t|* Am(t)) 
=|t|"*hm(t)(1 + o(1)) ast 30° 
<2] t[*-"E | ar | sex” (a — P). 
Now 
AV 
P(M< Vm & Piha | ay | *P | 1*-* dt 
0 
¥ 220-2e+2 Ean” [ap } =? 
a—P+] yor 
m 
22% -2p+2 5 x(a) Eee 
, 8 m n 
a—p+] yes 
But 
Safe). = _ 2a M2m-1 
or aL ; 
So 
irene eee ( -} ‘ex 
@—p)@—-P+H\d) ' 
kn \,_ 1 
: ( ae es 2me 2me 
( seers) s(1— =) 
e ) (log n) 
(by Lemma 2.2 and (4)) 
512 € kn \, l 
el ie 
oa (a—2mpo) (s- sme) 
( a, log n 
In 


(equation continued on p. 660) 
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660 
kn 
(". —> | for large n) 
tn 
e kn m l ——$—$—$$—$—————_— 
<171( +) ( m2 ec ee tome) 
(log n) . 
Ay | 
ee a | 
( ~ «-E 2 cat Os 
Choose Pp = lam for a fixed m and by this choice we have, 
P (Ty) < 171 ( i) oe ...(10) 
d } (logn)? * 


Then proceeding as in above we get 

r(«-p) a 
22%-2pt2 2 xm | ar : 
a—p+l prt 


m 


P(T2) < - 











But 
- hah: =<, | +- M2m-1 < 2M2m-1, 
So 
P (Ts) < 86 (4 ) Ee See 
4) ce ae (11) 
Now 
| Rm | = % ar Xr (w) x’ + a a’ X; (w) bi 
r | 
<|B ar X;y (w) iat ee ay Xr (w) an | 
Vn Vin 
S ata 2 = Vn 
Obviously, 


P(T) © P (Ty) + P (7) 


< 257 ea Ree 
d (log n)8 


(by (10) and (11)). 
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Then 
k 
£(G)= RACY, =k P(T) 
m= 


k 
< 257 5) ows 
ee: (4 (log n)3 


we (( ons) 


(by (4) and (8)). 
§ 3. Proof of the Theorem—Let 
Em = {w: Aom > Vom, Aom1 < — Vom4y} 
Fm = {@: Agm < — Vom, Aoma1 > Voms1} 
Dy, = {w: Aom > Vom, | Rom | < Vom} 
Dz = {o: Aomy1 < — Vomy1, | Rome1 | < Vomy1} 
D3 = {o: Agm < — Vom, | Rom | < Vom} 
D4 = {w: Aomy1 > Vomy1, | Romy | < Vomy1} 
and 
Ay = Dy (\\ Do, He = D3 1) Da. 


Obviously when H occurs, f (xem) > 0 and f (xom41) < 0, and also when He occurs, 
f (xem) < 0 and f (xom41) > 0. 
Therefore, if Hy Hz occurs, f (x) has aroot in (xam, xgm,1). Again we are to 
show that P (Em U Fm) > 0. 

Let Gm (x) and gm (t) be respectively the distribution function and the characte- 
ristic function of (Am/Vm). Then 

P (Em VU Fm) = (1 — Gm (1)) Geom, (—1) + Gom (1) (1 — Gom,1 (1). 
= (f3) 


Also 





gm (t) = exp | - EE ie x ar = h ee tlVm ) 
oe 1 


m 


= exp |- lr] Be | oft) ar x0 h(a, @|Vm ) 

x 1 

oe 
(equation continued on p. 662) 
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(1 + 0 (1) i (by lemma (2.1)) 
= exp {| —#| *-°@) 0°) (1 + 0 (1))} 
(by the definition of Vm). 
Therefore as m > ©°, gm(t) > exp (— | ¢ | *) in any bounded interval of t-values. 


Hence 


sup | Gm(x) — F(x) | = 0 (1). ..(14) 


So for « > 0, we have 


| Gam (— 1) — F(-— Li <8 


and 

| Goma (— 1) — F(-— YD) cae 
So 

Gom (— 1) > F(— 1) — «,  Gamai (— 1) > F(-—1)- « 
and 

1 — Gom (1) > 1 — F(l) -— «6, 1 — Gomy (1) > 1 — F(\)—«. 
Hence 

P (Em U Fm) > 2 (F(— 1) - e)\(1 —- F(1I)— 9 

> 2 F(— 1)() — F (1) as ¢ > O when m > © with n. 

Thus 


P (Em U Fm) = 35m = 8 > 0. . FS be 
Where 5 is an absolute constant. 
Let %m be the indicator function of the event Em U Fm. 
Then by (15) 

P {w: Ym = 1} = 8m and P {o: ym = 0} = 1 — Sm. 
Thus ym’s are independent random variables with 

E (nm) = 3m and V (ym) = (3m — 8° ) 


m<* 


It follows from Samal‘ (p. 439) that there are at least 


k ; 
, {k _ [>| ~ 2} pairs (Agp, Aop41) such that [ > 


+1qQp<spmr+ick. 
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If q be the number of such pairs, then 


q>3(e-[5]-2) 


< tk (for large k) 
log n 


S 32 tox (( 2) clog mor ) (by (4)). ..- (16) 


Let » = 2 ym, then for 0 < €; < 3m, by Chebysheff’s inequality 
@ 





V 1 
Pin-E@M|>qa< —% < + 
q? € Seti 


for large n, take 
1 


| log ( ( *_\ tog n)5!% )}rs 
flor (( dor ns pr 


i 





Then we have, 


P{|n-—E(M|Pqaty< 


log n 
(by (16)). 
Thus 
|» — E(») | <q «1 ona set S for which 
32 e {log (( re ) (log n)3/* ) fev 
P (Sy) < +S 2 17) 
Again y > E(n) —qa > q(@— 4). 
Taking 5 > 2 ¢ for large n, we have 
log n 
ee S. (9 oky Ne, 3) c, 
32 {log (( — ) (log n)5/* ) fsa 
tn 
log n 
ali (( ws re n) 5/2 ) 
og va g 
Hence 
i] ...(18) 


ee 
32 log ((—™ } dog m9 | 
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Now if S denotes the entire exceptional set 


Then 
P (S) & P(G) + P (Si) 


161 clog ((— ) dog m)>"* .. (19) 
\ tn 
Se eens 
logn 
Therefore from (18) and (19), we have 
9 peprogs pepeminbenmall FS 
32 log (( nt | (log n)5'* ) 
tn ] 


161 e log (( é. (log n)5!/* ) 
n 


< 
logn 
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GRONWALL, BIHARI AND LANGENHOP TYPE INEQUALITIES 
FOR DISCRETE PFAFFIAN EQUATION 


E. THANDAPANI 


Department of Mathematics, Madras University PG Centre, Salem 636011 
Tamil Nadu 


(Received 2 June 1987; after revision 18 July 1988) 


The object or this paper is to establish some new discrete inequalities of 
the Gronwall, Bihari and Langenhop type for Pfaffian equation. 


INTRODUCTION 


Itis well known that the discrete analogue of celebrated Gronwall’s lemma‘ 
and its generalization, Bihari’s lemma® established by Jones® and Sugiyama? (also see 
Agarwal and Thandapani!” and Pachpatte®) play an important role in the theorey of 
difference equations and numerical analysis. They have been used, for example, to 
obtain upper bounds of solutions as in Langenhop’ it is shown that lower bounds can 
be obtained similarly. The purpose of this note is to demonstrate how similar esti- 
mates can be derived for solutions of discrete Pfaffian equations and for the continuous 
case similar results are available in Grudo and Yarchuk®. 


Before giving the main result, we shall first introduce some notations which we 
t-1 

shall use throughout the paper. WN denotes the set {0, 1, ...}. The expression 2 6 (s) 
8=0 


represents a Solution of the linear difference equation A z (t) = b(t) forall t © N under 

the initial condition z(0) = 0...where A is the operator defined by Az(t) = 
“1 t-1 

z(t + 1) — z(t). It is supposed that 2 6 (s) = 0. The expression Lee (s) represents 
s=0 ag 

a solution of the linear difference equation z (¢ + 1) = c(t) z(t) forall t€ N with 


-1 
the initial condition z(0) = 1. It is supposed that Le c(s) = 1. We also define 
= 


A zr, (ti, fa) = Z(H + I, fe) — 2 (th, fa) 


A Zt» (ti, t2) = z (4, tg + 1) — z (hh, te) where the letters ‘1 and fg are 


used to denote the two independent variable which are members of N. 


MAIN RESULT 


Consider the discrete Pfaffian equation in two independent variables : 
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ty-1 ty -1 
x (ty, t2) = x (0,0) + =E fi (s1, to, x (si, t2)) + EF fo (0, 52, x (0, 52)) 
517° so 
yecta 
where 
fi (j = 1, 2) are real-valued and defined for all t1 > 0, te > 0 
and 
| f9 (t1, to, x) | < Fy (th, t2) Ws (| x1) (7 = 1, 2) sav) 


where the Fy are real-valued, non negative and defined for all 4% > 0, te > 0 and 
Wj (u) are continuous and non decreasing for u > 0 and Wj; (u) > Oforu > 0. 


Theorem —Let | x (ti, fe) | > O for t1 > 0 and te D> O. Then, if ty > 0, fe > 0, 
we have 


ty-1 to-1 
| x(t, te) |< fey DS Fy (51, t2) + af a(S Fe (0, s2)) |} 
8,=0 85-0 
=s-t3) 
where, 
Gow) = | BsO <u< a uo = | x (0,0) ] ) 4 
Wi (wa) , uo ; ...(4) 
“9 
and 
Eee (i = |, 2) are the inverses of Gi. 
Furthermore, 
eas i=1 
Ix > 61 Sri, 9) +a] o(- > n@ s2)) |} 
8,=0 is 85-0 
vat hs 
for all ty > 0, to > O for which 
t,-1 to-1 
x a 
Ss Fy (51, t2) + Gy |G; f= SS Fe (0, s2)) | € dom (G;") 
8,70 8-0 


and 
to-1 


— > Re, 5) € dom (G;}). 


85~0 


DISCRETE PFAFFIAN EQUATION 667 


Proor : By putting | x (t1, f2) | = u (ty, t2), we conclude from (2) that 
Fs t-1 | Pea | 

u(t,t2) <u~+ ; Fy (1, t2) Wi (u (si, t2)) + Z Fe (0, 52) 

7 =0 


#9 
We (u (0, s2)). ...(6) 

Now put 

t,-1 ty-1 

R (ti, t2) = uo + ZB Fy (51, to) Wy (u (51, t2)) + 2 Fe (0, sa) 
8,=0 85-0 
Ws (u (0, s2)) 

to obtain 


AR = Fi (th, ta) Wi (u (hy, f2)) 


< Fy (h, te) Wy (R (h, t2)) 


since W, (u) is non decreasing and wu (ft, t2) & R(t, tg). Since R(t, t2) > 0 and 
Wy (u) > Oforu > 0 foru > 0, this inequality together with (4) yields 
ty-1 
Gy (R(t, t2)) — Gi (R (0, f2)) < a Fy (51, f2). -.(7) 
re 
1 


From (6), we have, on putting 11 = 0 
to~1 
u (0, t2) {uo + tae Fy (0, s2) We (u (0, s2)). 
2 
As before, we conclude that 
t2-1 
Go (R (0, t2)) — Ge (R (0, 9)) < -s F2 (0, 2). 
2 


Since R (0,0) = uo and Gz (uo) = 9, we have 


to-1 
R02) < 6; [ > O50). (8) 
85-0 
Hence, since Gi (R(t, f2)) > Gi (u(t, to)), we conclude from (7) and (8) that 
t.-1 | 


2 1 
Gy (u (t1, t2)) < Gi 19," DS Fo (0, =) | + s: Fy (si, t2) 
85-0 8,~0 


and the result (3) follows. 


We now prove (5). Again using the notation | x (t1, t2) | = u(t, t2), We find from 


(1) and (2) 
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at 
u(t, t2) > u(s1, 52) — 2 Fi(n, te) Wi (u(r, to) 
Naa 
tol 
— ZL Fe(si, ro) We (u (si, 72)) aca 
To 89 


fort) > 0, t2 >OandO xcs GUKi= ay 


Let 
ty-1 
R (si, $2) =u (t,t) + D Fi (ri, te) Wi (u(r, 2) 
oh Bs 
to-1 
+ LZ Fe (si, r2) We lu (1, 72). 
To~89 
Then 
ARs, = — Fo (s1, 52) We (u (s1, $2) 


and since R (51, s2) > u (si, $2) 


ARs, > — Fo (s1, 52) We (R (51, $2)). 


Thus 
-1 
Ce (R (s1, t2)) — Go (R (51, 52)) > ar Ba Fo (s1, re) 
2789 
and we conclude that 
ty-1 
Ge (R (s1, t2)) > Ge (u(s1, 52)) — 2B Fo (51, 12). 


To % 
This implies that 
ty-1 
R (s1, t2) > G5'[Ge (wu s1, $2)) — SS) Fe (s1, r2)]. ...(10) 
"989 


Putting so = f2 in (9), we obtain 
ty-1 
u(t, to) > u(s1, f2) — 2 Fy (ri, te) Wi (u (rn, t2)). 
fe i | 


Hence 


t,-1 
G1 (R (t, f2)) — Gy (R (51, t2)) > — 53 Fy (ri, te). 
of Be | 
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Since R (%, f2) = u (ty, te), (10) implies that 
t-1 
Gy (u (41, t2)) > > Fy (ry, t2) + Gy | G; | Gs (u (s1, 52)) 
t.~1 


— 2 Fe (51, re) } 
To" 85 


and putting s; = 0, so = 0, we obtain 


‘71 ty-1 


Gx (u(t) > — SF (1, ) + G 1% | - > Fx (0, 72) |} 


1,70 
This proves (5). 
Remark 1: If W; (u) = We (u), (3) and (5) simplify : 

to-1 

| x (41, f2) | < Ga > Fy (51, t2) + SS Fe (0, se) ] 
8570 
id 

x(t, t) | > G7 (5 Fates) = Ss Fs (0, 2) J. 
85=0 


Remarks 2: If Wi (u) = We in = u, we have 


Pid | t,-1 
| x (41, 2) | & | x (0,0)| TW [1 + Fi (81, t2)) ip [1 + Fe (0, s2)) 
8,0 97 
t- 1 t,-1 
<|x(0,0)|exp[ = Fy, (si t2) + 2% Fe (0, 52)}. 
8,70 §5=0 
t,-1 t.-1 


| x (ti, t2)| > | x(0,0)| HM [1 — Fy (81, t2)) MI [1 — Fe (0, sa) 
: 8,=0 


= t 71 
aT 1 


> |x(0,0)|exp[— = Fi(si, te) - & Feo (0, 52)). 
8,0 8570 
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Hans H. Storrer® has introduced the notion of critical left ideals in rings and 
has given a generalisation of primary decomposition to non-commutative 
rings. In this paper critcal left ideals in Near-rings are introduced and the 
relation between critical left ideals and left ideals of type O is studied. 
Further primary decomposition is also exhibited in modules over neuir-ring 
under some conditions. 


1. INTRODUCTION 


All near-rings are assumed to be zero-symmetric right near-rings with identity. 
Throughout this paper near-ring under consideration is denoted by K. The definitions 
of K-module, K-subgroups and submodules are as given in Pilz?. However for the sake 
of continuity the definitions are given again. 


Definition 1.1—Let (M, +) bea group and K be a near-ring such that there 
exists a mapping u: K x M — M satisfying the conditions. 


(kK+k')m=km+k'm 
(kk') m = k (k' m). 


l.m = m for all k, k’' © K,m € M and | is the identify of K. Then (M, +, u) 
is called a K-module. 


Definition 1.2—A subgroup N of K-module M is said to be a K-subgroup of @ 
if (N, +) is a subgroup with KN C N. 


Definition 1.3—A normal subgroup N of M is called a submodule of M if 
k(m+n)—kmeEN. Forallm€ M,n€ Nandk E K. 


The concepts of essential extension and rational extension were introduced by 
Barua! for near-rings. Here we give a variant of the above notions similar to modules 
over a ring. 


Definition 1.4—Let M be a module over a near-ring K. M is said to be (a) a 
‘module essential extension’ of a non-zero K-subgroup N if for every non-zero sub- 
module N’ of M, NQ N’ < 0; (b) an ‘essential extension’ of a non-zero K-subgroup 
N if for every non-zero k-subgroup N’, N (1 N’ # 0. 
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Clearly essential extension implies module essential extension. The converse need 
not be true. 


Example—Let K = {0,a,b,c} be the klien group under addition and the 
multiplication be defined by the following table. 





0 a b c 
0 0 0 0 0 
0 a 0 a 
0 0 b b 
c 0 a b c 


owe ce 





SS ee RS A 





Then M = K is module essential extension of N = {0, a} but not an essential exten- 
sion of N. 


Definition 1.5 : Rational Extension—Let N be a non-zero K-subgroup of a K- 
module M. M is said to bea rational extension of N if the following condition is 
satisfied : 


NCA C M, Aisa K-subgroup of M and f: A > Misa 
K-homomorphism such that ker f D N implies f = 0. 


In modules over rings it is known that a rational extension implies essential 


extension. This need not be the case in case of near-ring modules as seen from the 
following example. 


Example: The dihedral group Dg = {0, a, 2a, 3a, b,a + b, 2a + b, 3a + b}. 
With the addition and multiplication operations defined below is a near-ring with 


identity a. It is a rational extension of each of its non-zero K-subgroups but not an 
essential extension of its non-zero K-subgroups. 


ae 0 a 2a 3a b a+b 2a+b 3a+b 
0 0 a 2a 3a b a+b 2a+b 3a+b 
a a 2a 3a 0 a+b a+b 3a+b b 

2a 2a 3a 0 a 2a+b 3a+b b a+b 
3a 3a 0 a 2a 3a+b b a+b 2a+b 
b b 3a+b 2a+b a+b 0 3a 2a a 
a+b a+b b 3a+b a+b a 0 3a 2a 
2a+b a+b ath 4 3a+b 2a a 0 3a 
3a+b 3a+b 2a+b atb b 3a 2a a 0 
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0 a 2a 3a b a+b 2a+b 3a+b 
0 0 0 0 0 0 0 0 0 
a 0 a 2a 3a b a+b 2a+b 3a+b 
2a 0 2a 0 2a 0 0 0 0 
3a 0 3a 2a a b a+b 2b+b 3a+-b 
b 0 b 2a b b a+b 2a+b 3a+b 
a+b 0 a+b 0 3a+b 60 0 0 0 
2a+b 0 2a+b 2a 2a+b b a+b 2a+b 3a+b 
3a+b 60 3a+b 0 a+b 0 0 0 0 





Definition 1.6—A K-module M is said to be (i) ‘Uniform’ if it is an essential 
extension of each of its non-zero K-subgroups; and 


(ii) ‘Strongly uniform’ if it is uniform and is a rational extension of each of its 


non-zero K-subgroups. 


Notation: If X is a subset of a K-module M then <X> stands for the sub- 
module of M generated by X. 


We assume that K-module M satisfies the property (p): 
“=< Ny (\ No > = < Ny > N < Ne > for any two K-subgroups Ny and Ne of M”’. 


Any near-ring K in which every K-subgroup is a submodule of K satisfies this 
property. 

Proposition 1.7—Let M be a K-module with ascending chain condition (ACC) 
on K-subgroups and satisfying the property (p). Then M has a submodule which is 
uniform. 


Proor : Suppose M has no submodule which is uniform. Then there exists non- 
zero K-subgroups Ni and N| such that Ny ON, = 0. Therefore M1 1M, = 0 


where My = < Ny >and M, = <N, >. By hypothesis Mj, is not uniform. So 
there exists non-zero K-subgroups No and N, of M, such that No 1) N, = (0). There- 


fore M2 1 M, = () where < No > = Meoand < N, >= M, . Repeating the 


argument, we get a sequence {Mn} of submodules of M which are not uniform such 
that Mi © My + M2 C Mi + Mo + Mg C ..contradicting the A.C.C. on K- 
subgroups. 


Hence M has a submodule which is uniform. 


2. CRITICAL Lert IDEALS 


Following the notion of critical left ideal as given by Storrer? the concept of 
critical left ideals in near-rings is introduced. 


674 K. YUGANDHAR ET AL. 


Definition 2.1—A left ideal P of a near-ring K is called a critical left ideal if 
M = K/P is a strongly uniform K-module. 


Definition 2.2—(a) A left ideal L of K is of type 2 if M = K/L is a K-module 
which has no nontrivial K-subgroups. 


(b) A left ideal L of Kis of type Oif M = K/L is a K-module which has no 
nontrivial submodules. 


From the definition it follows that every left ideal of type 2 is a critical left 
ideal. However we can still strengther the result as follows. 


Proposition 2.3— Every left ideal of type O of K is a critical left ideal. 


Proor: Let L be any left ideal of type O of K. Suppose At = Li /L, i = 1, 2 
are two non-zero K-subgroups of M = K/L. Assume that Zi € Lo and Lo € L1. 
Choose r; € Li\ le andre € Lo\ Ly. Since LZ is a modular left ideal of type O, 
(K: L) = P is a O-promitive ideal of K and M = K/L isa K /P-module of type O. 


By density theorem (Pilz”, 115) there exists n + p€ N/P such that (n + p) 
(ro+ L)=1r1+ L which implies that mro—r1 © L C Ly. Consequently 
nro © Li (1) Lo and nrg & L, Thus nrg + LE Ay A Ag. This shows that inter- 


section of any two non-zero K-subgroups of K/L is non-zero. Thatis M = K/L is 
uniform. 


Now to establish that 4 = M/L is a rational extension of each of its non- 
zero K-subgroups, assume that there exists two non-zero K-subgroups A; and Ao of M 
such that A; C A» anda K-homomorphism f: Ag > M with ker f D Aj. As before 
suppose Ai = Li/L,i= 1,2. If Ly = Le, then clearly f = 0, Otherwise choose 
any re © L2\ Lj; andr; € LiNL. Considering M = K/L as a K“P-module and 
appealing to density theorem again, there exists n + p © N//P such that (n + p) 
(ri + L) = ro + L which implies nry — ro € L. Since fis a K-homomorphism and 
nry + L & 41 we have f (re + L) = f{nry + L) = 0. 


This implies that f = 0 and M ina rational extension of Aj. Hence Z is a cri- 
tical left ideal of K. , 


In the class of all critical left ideals of a near-ring K, we introduce a relation as 
follows. 


Definition 2.4—If 1 and J are two left ideals of K, / is said to be related to J if 
there exists a © J, b © J such that 
Ja“) = Jb-1 where Ja-1 = {r € K /ra € /}and 
Jal= {re K/rb € J}. 


Lemma 2.5—Let K be a near-ring and P a 


nd Q be two critical left ideals of K 
then the following statements are equivalent : 
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(i) P is related to Q. 


(ii) A non-zero K-subgroup of K /P is isomorphic to a non-zero K-subgroup of 


K/Q. 
Proor : Suppose P and Q are related, then there exists a and b such that 
Pal= Qblanda & P,b & Q. 
Then Pa-} is a left ideal of K and K /Pa-1 is a non-zero K-module Define : 


b:K/Pa1—>K/P 
by 
v(x + Pa!) = xa + P. 


The map ¥ is an isomorphism of K// Pa-! onto a K-subgroup L/P of K/P. Similarly 
K/Qb~1 is isomorphic to a K-subgroup L’,/Q of K,/Q showing that L “P ~ L'/Q 
if Pa-1 = Qb-1, 


Conversely suppose a non-zero K-subgroup NP is isomorphic to a non-zero 
K-subgroup M/Q of K /@. 


Let ¢ be a K-isomorphism from NV /P to MQ. 
Since ¢ = 0, there exists d (a + P)=b+O0 0. 
¢ [r (a + p)] = $ {ra + P) 

r¢ (a+ P) 
=r(b+Q)=rb+Q. 


Since ¢ is mono, ra € P <= rb € Q that is Pa-1 = Qb-1, 


I 


Proposition 2.6—The relation P is related to Q is an equivalence relation in the 
class of all critical left ideals. 


Proor: The only condition to be verified is the following: If /, J, Z are cri- 
tical left ideals such that / is related to J and J is related to L then / is related to L. 


By Lemma 2.4, there exists non-zero K-subgroups My, /J, Ni /J of K/J and 
K/J respectively such that M;,/J = N;/J. Similarly there exists non-zero K-sub- 
groups No /J and M2 /L of KJ and K ’ L respectively such that No/J = Me/L. 
Since K/J is uniform, we have Ni M Neo = wf 


Put N = Ni © Ne. Then N/J is isomorphic to a non-zero K-subgroups of 
K /I and also NJ is isomorphic to a non-zero K-subgroup of K/L. Hence a non- 
zero K-subgroup of K//J is isomorphic to non-zero K-subgroup of K/L. Hence /is 


related to L. 
The equivalence class containing P is denoted by [P]. 
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3. ASSOCIATED LEFT IDEALS OF A MODULE M 
Definition 3.1—A_ critical left ideal P is said to belong to M if there exists 
0 # x € M such that Ann (x) = P. 


Here we identify, the critical] left ideals related to P and say [P] is associated to 
M. With this identification, the set of all critical ideals belonging to M is denoted by 


Ass M. 


Theorem 3.2—If M satisfies A.C.C. on K-subgroups and also satisfies property 
(p). Then there exists a non-zero K-subgroup B of M which is strongly uniform. 


Proor : By Proposition 1.7, there exists a non-zero K-subgroup N of M which 
is uniform. 


If N is strongly uniform, nothing to prove. Suppose now NJ is not strongly uni- 
form. But by our choice AN is uniform. 


Therefore N cannot be a rational extension of each of its non-zero K-subgroups. 


That is, there exists a maximal K-subgroup A of Nsuch that Nis nota rational 
extension of A (since ™ satisfies A.C.C. on K-subgroups). Hence there exists K-sub- 


group N’ such that A C N’ C N,and¢: N’ + N,a non-zero homomorphism with 
= 


Ker ¢ D A. 
Let ¢ (V’) = B+ 0. 


Then B is a non-zero K-subgroup of N. We claim that B is a rational extension of 
each of its non-zero K-subgroups. 


Let B’ 40 bea K-subgroup of B. 
Consider 4-1 (B’) = A’ then MS A’ C N and Nis a rational extension of A’ (by the 
maximality of A). 
Let f: B° + Bbea homomorphism where B’ C B" C Bwith Kerf 2.3’. 

If $1 (B’) = A” then ¢-1(B) > ¢-1 (3") "SD FBS. 
Then f o ¢: A” >N isa homomorphism such that Ker f°¢DA’. But Nisa 
rational extension of 4’. Therefore tf °¢=O0and hence f = o on B”. 
That is, B is a rational extension of B’, 


Therefore BC Nand B isrational extension of each of its K-subgroup. As JN is uni- 
form, B is also uniform. Therefore B is Strongly uniform. 


AS a corrollory to the above, we have : 


Corollary 3.3—If M sati x 
(pers atisfies AC.C. on K subgroups and also property (p), then 
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PRroor: By Theorem 3.2, Mhasa K-subgroup N which is strongly uniform. Let 
0 ~x EN. Then Kx = K/YP where P = Ann (x) and Kx is strongly uniform. 
Hence P is a critical left ideal associated to M and Ass (M) ~. Some of the proper- 
ties are the following : 


Proposition 3.4—Let M bea K-module satisfying A.C.C. on K-subgroups. 
(a) If M is the union of submodules of Mi, then Ass M = U; Ass Mi, 
(b) If P is a critical left ideal then Ass (K/P) = {[P}}. 
(c) If N C M, then Ass N C Ass M C Ass N \) Ass M JN. 
(d) If M is the direct sum of submodules Mi, then Ass M = Ui Ass Mi 
Proor : (b) To show Ass ‘K_/P) = {[P]}. 
Consider Ann (1 + P) =[r © K/r(1 + P) = 0) = P. 
If there exists x € K such that Ann ear) == 0, 
Q.1-1 = Ann(x + P) =[re K/Yr (x + P) = 0] = Px-l, 
Hence P is related to Q. 
Therefore Ass (K./P) = {[P]}. 
(c) If NC M, then Ass (N) C Ass (M) is clear. 
Let N be a submodule of M. 
Let P € Ass M and P = Ann (x) Oseex CM. 
Case (i)—If Kx \. N = <0> 
Ann (x + N)in M/N = [ye K/y(x + N)=0 + NJ] 


=[y € K/yx = 0] 
= Ann (x) in M. 


That is P € Ass(M Y/N). 
Case (ii)—If Kx 1. N+ <0>. Let N’ = Kx (1) N consider N’ as K-module. 


Then Ass (N’) # © by corollory 3.3. 
Let Q € Ass(N’) C Ass N. 
QO € Ass (N’) > There exists y € N’ such that Ann (y) = Q 
y € Kx =» There exists a € K such that y = ax. 
Q=Ann(y) = {z/zy = 0} 
= {z/(za) x = 0} 
{z/za € P} = Pa“l, 
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So Pa! = Q1-! and P is related to Q. 
Hence [P] € Ass (N). 


4. PRIMARY SUBMODULES 


A K-module M is said to be coprimary if Ass M consists of a single clement and 
and a submodule N of M is said to be primary if quotient module M/W is co- 
primary. 


Theorem 4.1\—Suppose 4/ is a module which satisfies A.C.C. on K-subgroups 
and uniform. Then M is coprimary. 


Proor: Let P € Ass M, P = Ann (x), Kx = K/P. 
If OQ € Ass M, then Q = Ann(y), Ky = K/Q. 
Since M is uniform, Kx } Ky + 0. 


Let C = Kx 1 Ky. Then C is isomorphic to a non-zero K-subgroup of K_“P and also 
a non-zero K-subgroup of K/ Q. 


Therefore P is related to Q. 
O~ PE Ass M. 
Therefore Ass M = {[P]}. 


Definition 4.2—Let N be a submodule of M. N is said to be irreducible if for 
any two submodules N; and No, N Cc Ni,NC No>WN & Ni 0) No. 
7 = 


And N is said to be ‘strongly irreducible’ if for any two K-subgroups Ny and No, 
NC NandNC Ne >+NC NA No. 
a a of: 


If M satisfies property (P), clearly a submodule JN is irreducible if and only if it 
is strongly irreducible. 


Hence if JN is an irreducible submodule then M “N is uniform. 


Theorem 4.3—If M is a K-module which satisfies A.C.C. on K-subgroups then 
every submodule of M is a finite intersection of irreducible submodules. 


This can be proved by using the A.C.C. on K-subgroups. 


Definition 4.4—A primary decomposition of a submodule N of M is a represen- 
tation of N as a finite intersection of primary submodules. 


Theorem 4.5—If M is a K-module satisfying A.C.C. on K- 


subgroups and pro- 
perty (p), then every submodule N of M has a primary decomposi 


tion. 
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Proor: By Theorem 4.3, a submodule N of M can be written as N = Ny -) N2 


1) ... (\ Ne where each Ni is irreducible submodule. By the property (p) each Ni is 
Strongly irreducible. 


Then for each i, M,N; is uniform and satisfies A.C.C. on K-subgroups and it 
is coprimary by Theorem 4.1. Thus, Né is primary. 


5. Z-S-CopRIMARY SUBMODULES 


It can be shown that every primary decomposition has reduced primary decom- 
position and uniqueness by the familiar methods. 


Lemma 5.1—Let K bea near-ring with | satisfying the conditions. 
(a) Possesses A.C.C. on ideals. 
(b) Every left K-subgroup is an ideal of K. 
Then every critical left ideal of K is a prime ideals. 


Proor : By condition (b), P isan ideal of K and by definition of critical left 
ideal, K / P is strongly uniform. 


If P is not prime, there exists ideals A and B such that P = Ar. 3 B and 
BA C P. Since P A, there exists a € A anda a ole age 
Define mapping f: K.P > KP as follows : 
i (xet+ £).=9xa + P 
clearly this is a module homomorphism and Ker (2.87 -P and f= 0. 
For if f = 0, then xa © P forevery x € K. 
That is Ka C P which is not the case since a € Ka anda & P. 


This is a contradiction to the hypothesis that KP is a rational extension of 
each of its non-zero K-subgroups. 


Hence P is prime. 

Definition 5.2—An ideal Q of K is said to be Z-S-coprimary if Px CQ>PCr(Q) 
orxE Q 
where r (Q) is the intersection of prime ideals. 


Theorem 5.3—Ilf K-satisfies conditions above, then an ideal Q which is Z-S- 
coprimary is coprimary in the sense of this article. 


Proor: P € Ass(K “Q) and P = Ann (x + Q). 
That is Px © O,x & Q. 
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Therefore P C r(Q) but r(Q) C P (since P is prime) 
therefore P = r (Q) which implies Ass (K/Q) = {[P]}. 
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_In the present paper I have obtained some of the properties of an almost 
semi-invariant submanifold of an §P-Sasakian manifold. The intergrability 


conditions of the distributions D, D+, D, D @ {#}, Di @ {€} and D @ {8} 


have also been discussed. 


1. INTRODUCTION 


~~ 


Let M be ann dimensional C®-manifold. If there exists in atensor field F 
of type (1, 1) a vector field € and a 1-form » satisfying 


F2 X¥ = X — y (X) &, 1 (E) = 1. (iT 


Then M is called an almost paracontact manifold. 


Let g be a Riemannian metric satisfying 


n (X) = g (X, &) me dey)’ 
n (FX) = 0, FE = 0, rank (F) =n — | etlss) 
g (FX, FY) = g (X,Y) — n (X) 1(Y). (1.4) 


Then the set (F, &, », g) is called an almost paracontact Riemannian structure and 
the manifold is called an almost paracontact Riemannian manifold2. 


Moreover if we define 


'F (X, Y) = g (FX; Y) met 6 BEY | 
then 

eae poor (Y.- x) ...(1.6a) 

'F (FX, FY) = 'F (X, Y). ...(1.6b) 
Now, let us consider a manifold in which the |-form 7% satisfies 

(vx 1) (Y) — (vy 1) (X) = 0 (1.7) 


(Vx Vr 1)(Z) = (—8 (X, Z) + 1(X) 9 (Z)) 1 (Y%) 
+ (—g (X,Y) + n(X) 1 (¥)) 1(Z) rie) 
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where VY denotes the covariant differentiation with respect to g. Furthermore, if we 
put 

n(X) = 8 (X, &), vx & = FX (1.9) 


then it is easily verified that the manifold in consideration becomes an almost para- 
contact manifold. Such a manifold is called a p-Sasakian manifold. 


If the 1-form y in M satisfies 
(vx n)(Y) = — (X,Y) + n(X) 4 (¥) ...(1,10) 


we can easily show by putting y (Y) = g (X, &) and (vx y)(Y) = ‘F(X, Y), that the 
manifold is p-Sasakian. Such a manifold is called an sp-Sasaksan manifold’. Thus 
in an sp-Sasakian manifold, we have 


F(Y,Z) = — g(Y, Z) + » (Y) »(Z) Reet gh Be 
and 
(vx 'F)(Y, Z) = (— g (X,Y) + 4 (X) n (Y)) n(Z) 
+ (— g (X,Z) + 1(X) n (Z)) 0 (¥). (1.12) 
Let M be an m-dimensional submanifold immersed in an sp-Sasakian manifold 


M. Let TM and TM! be respectively the tangent and the normal bundle to M. 
Suppose the structure vector field & is tangent to M and denote by {&} the 1-dimensional 
distribution spanned by — on M and {}4 the complementary orthogonal distribution 
to {§} in TM. Foreach ¥ € [ (TM), put 


FX= bY + c¥ (1.13) 
where bX € [ ({E}1) andcX¥ ET (TM)+. Thus dis an endomorphism of the tangent 


bundle 7M and c is a normal bundle valued |-form on M. 


Definition 1.1—The submanifold M of the sp-Sasakian manifold M is said to be 
an almost semi-invariant submanifold if its tangent bundle 7M has the decomposition 


TM =D®D1@D@ {3} .. (1.14) 
where 
(1) Dis invariant distribution on M, i.e. 
F (Dz) = Dz 


(2) D1 isan anti-invariant distribution on M, i.e. 


F (D2) Cet et forx € M. 


SEMI-INVARIANT SUBMANIFOLD 683 


(3) D is neither invariant nor an anti-invariant distribution on M, i.e. b Xz #0 
and C ¥z; = 0 for anyx € M and Xz € Dz. 


(4) {&} is the distribution spanned in M by the vector field &. 


2. Basic RESULTS 


Let M be an almost semi-invariant submanifold of an sp-Sasakian manifold M. 
We denote the Riemannian metrics of M and M both by g. Let P,Q and L be the 


projection morphisms of TM to the distributions D, D+ and D respectively. Then for 
X € T(M), we have 


X = PX+QX+4+ LX + n(X)E. beh 


Now, we take YET (D). Then bX¥ 40, CX¥ 0. Thus C defines a vector sub- 
bundle CD :x—> CD, of TM:. 
For any N€ I'(TM!), we put 
FN =tN+{fN (202) 
where tN and /N are respectively the tangential and the normal components of FN. 
Then we have, 


Eir De. CD) = 0. sexleea) 


Next, we denote by v the orthogonal complementary vector bundle to FD4 ®© CDin 
TM-+. By (1.4), we have 


BPX CY) mig as FI) =. 2(X, 1) = 0 MX ET (D+) 


Ye YI (D). ...(2.4) 
Thus 
TM1= FD+ @ CD @ v. es oie) 
Lemma 2.\—The morphisms ¢ and f satisfy 

t(TM+) = D* ® D ...(2.6) 

t(FD!) = D* » (2-4) 

D) =D (2.8 

t(C D)=D (2.8) 

D D (2.9 

f (CD) = CD. (2.9) 
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Proor : Let N€ I (TM+), then 
g (tN, X) = g (FN, X) = g (N, FX) = 04 X ET (D) 
and 
g (tN, 5) = 8 (FN, &) = g (N, Fe) = 0. 
Thus ‘N € I (D+ @ D) and we get (2.6). Next for each ¥ € I (D*), we have 
X = F2X = tFX + CFX = tFX 
which implies (2.7). The proof of (2.8) and (2.9) is same as given in Bejancu and 
Papaghiuc!. 
Definition 2.1—An almost semi-invariant submanifold M in an sp-Sasakian 
manifold M is said to be a semi-invariant submanifold if we have D = {0}. 
Proposition 2.1—Let M be an almost semi-invariant submanifold of the 


sp-Sasakian manifold M. Then the endomorphism b: TM > TM is a para-/ structure 
on M, that is, 63 — b = 0 if and only if M is a semi-invariant submanifold. 


ProoF : From (1.13), we see that 


(65 — b)X =0 forany YE r(D @ D+ @ {&}). 


Since b Dr = Dz, we see that 


(b8 — b) (D) = {0} 
if and only if 
(62 — 1) (D) = {0}. 
Applying F to (1.13) and using (1.13), (2.2) and (1.1) we get 
X= F2X¥ =X + ChY¥Y + tCX + fCX for X € T (D). 
Equating the tangent part, we get 
(62 — 1) = — tC. 
Therefore, 
tC (D) = {0} 
which with the help of (2.8) gives D = {0}. Hence the submanifold is semi-invariant. 
Proposition 2. 2—Let M be an almost semi-invariant sub-manifold of the 
sp-Sasakian manifold M. Then M is a semi-invariant submanifold if and only if 


‘ge ik tl (2.10) 
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Proor: We see that if NE (FD+), we have fN = 0 and for NE I (v), 
JN = FN. By (2.9) wesee that fis an automorphism on CD. Hence (f* — f) 


(CD) = {0} if and only if 


6 seit —ae 9 (CD) = {0}. male 
By means of (1.13), (2.2) and (1.1), we get 
CX = FA(CKX) = b/CY + CY'+ C + Cx + f2CX pod e2l 2) 


for any Y € (CD). 
Thus (2.11) is equivalent to 

Ct (CD) = {0} ral2. i) 
which with the help of (2.8) gives (2.16). 


Lemma 2.2—Let M be an almost semi-invariant submanifold of an Sp-Sasakian 


manifold M. Then we have 


(b2 + tC) X = X — n(X)E, (Ch + fC) X¥=0 (2.14) 
(bt + tf) N =0, (f2 —1+ Ct)N =0 (2.15) 
Ue f 4 Cif) iN 0 (2.16) 
(b3 — b + 1Cb) X = 0 (2.17) 


for any X € [' (TM) and N € ['(TM?). 


Proor : It follows directly from (1.1), (1.13) and (2.2). 


Let vy (resp Y) be the Riemannian connection on M (resp M) with respect to the 


Riemannian metric g. The linear connection induced by Vy on the normal bundle 7M+* 
is denoted by y*, then the equations of Gauss and Weingarten are respectively given 
by 


VeY = we + h(X,Y) (2.18) 
and 
vxwNe= —AnX+,N (2.19) 


for all ¥, Y € [.' (TM) and N € ['(7M*), A is the second fundamental form of M 
and Aw is the fundamental tensor with respect to the normal section N and 


g(h (X,Y), N) = 2 (Av X,Y). _. (2.20) 
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Lemma 2.3—Let M be an almost semi-invariant submanifold of the sp-Sasakian 


manifold M. Then we have 


P (u (X, Y)) = FP (yx Y) — 4 (Y) PX etzeel) 
Q (u (X, Y)) = Ot (A(X, Y)) — 1 (Y) Ox ey ee 
L (u (X, Y)) = bL (yx Y) + Lt(h(X, Y)) — 1 (Y) LX 022.23) 
n(u (X, Y)) = — g (FX, FY) .-.(2.24) 


h (X, FPY) + h(X, bLY) + vy. FOY + yy, CLY 


= FOYx Y+CLyx¥+/(h(X, Y)) BA ees) 
where 
u(X,¥) = yx FPY + yx bLY — Aroy X — Acry X wheeee) 
for all X¥, Y € T (TM). 
The proof follows from Bejancu and Papaghiuc!. 


Lemma 2.4—Let M be an almost semi-invariant submanifold of an sp-Sasakian 
manifold M. Then we have 


vx § = FX, h(X,§) = 0; for any X¥ € I'(D) 0 (2:27) 
Vr §=0, h(Y, &) = FY: for any Y € ['(D+) ...(2.28) 
vz& = bZ, A(Z,8) — CZ, forany Z € T(d) (2,29) 
ve& = 0, h(E, &) = 0. .. (2.30) 


PrRoor : We have 


VrE= rE + A(X, 5) (2.31) 
which with the help of (1.9) and (2.1) gives 

Vx & + h(X,&) = FPXY + FOX + bLX + CLX vae( aioe) 
forany YET (TM). 
Now, (2.27)-(2.30) follows from (2.32), 


Lemma 2.5—Let M bean almost semi-invariant submanifold of the sp-Sasakian 
manifold M, then we have 


Arx Y + Ary X = 0. 
for all X, YE lr (D+), 


---(2.33) 
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Proor : With the help of (1.4), (2.18), (2.20), we have 
g (Arx Y, Z) = g (h(Y, Z), FX) = g (92 Y, FX) = g (Fz, X) 


= § (yz FY, X) = — g (FY, vz X) = — g(h(X, Z), FY) 
= — g (Ary X, Z) 
forall X,Y E [(D*‘)andZET (TM), which implies (2.33). 
Lemma 2.6—Let M be an almost semi-invariant submanifold of the sp-Sasakian 


manifold M. Then we have 


veVET(D) forany VET (D) ... (2.34) 
veVY EF (D") forany V € F(D+) (2.35) 
veWeET(D) forany W € TD). orl 2zo0) 


The proof follows from Bejancu and Papaghiuc!, 
Corollary 2.1—Let M be an almost semi-invariant submanifold of the sp-Sasakian 


manifold M. Then we have 


[X,&] E€ [(D) forany ¥ € [(D) et230) 
[Y,§)E€ [ (D+) forany YET (D3 »--(2.38) 
[Z, J € r (D) forany Z€ r (D) ...(2.39) 


The proof follows from Lemmas (2.4) and (2.6). 


3. INTEGRABILITY OF DisTRIBUTIONS 


Theorem 3.1—Let M be an almost semi-invariant submanifold of an sp-Sasakian 
manifold M. Then the distribution D is integrable if and only if 
h (X, FY) = h(Y, FX). Fg 
ProoFr : Wy using (2.27), we have 
g ([X, Y], &) = g (vx Y — vy X, &) 
m — (9x 5, Y)-+...4- 9 (vr 6, X) 
= — g (FX, Y) + g (FY, X) 
= 0 


for all X, Y € [ (D). 
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Next, from (2.25), we have 
h (X, FY) = FO vx Y+ CL yx ¥ +f (h(X, Y)) reoael 
for any X,Y € ['(D). Hence we have further 
h (X, FY) — h(Y, FX) = FQ ([X, Y]) + CL (LX, Y) ey 
which proves the theorem. 
From Theorem 3.1! and (2.37) follows : 


Corollary 3.1—The distribution D @ {&} is integrable if and only if (3.1) is 
satisfied. 


Theorem 3,2—The distribution D+ is never integrable. 
ProoF: For X¥,¥ € ['(D+), (2.26) gives 
u (X,Y) = — Ary X. 
Applying F to (2.21) and using (1.1), we get 
P(vx Y) = FP (Ary X), forany X,Y € [(D+) 
which with the help of Lemma (2.5) gives 
P ([X, Y]) = FP (Ary X — Arx Y) = 2FP (Ary X) 
showing the non-integrability of D1. 
From Theorem 3.2 and (2.38) follows : 


Corollary 3.2—The distribution D1 @ {€} is never integrable. 


Theorem 3.3—The distribution D is integrable if and only if 


Acy X — Acx Y — yx bY + yr bX ET (Die D ® {5}) ...(3.4) 
and 


h (bX, Y) — h(X, bY) + Vy CX — vicrer (CD + v) ---(3.5) 
forany Y, YET (D), 


Proor: ForanyX,YEPr (D), using (2.29), we get 


#0, ¥10) = 8 (vx ¥ — wr ¥,8) = (X, v8) — g(V, yr 8) 


= & (X, bY) — g (¥, bX) 


=a. ...(3.6) 
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Now, for any Y, Y € [ (D), (2.26) gives 


u (X,Y) = Vx bY — Acy X. ee ED 
Applying F to (2.21) and using (1.1) and (3.7), we get 

P (Vx Y) = FP (Acy X — yx bY) SAAD, 

P ([X, Y]) = FP (Acy X — Acx Y — yx bY + yy bX) = 0 .-.(3.9) 


if and only if (3.4) is satisfied. Applying F to (2.25) and taking the components in 
D1, we get 


QO vx ¥ = Of (h (X, bY) + yi, CY — fh(X, Y)) 
which further yields 


Q (LX, Y]) = Qt (A(X, bY) — h(Y, bX) + yi. CY — yb CX). 


Hence D is integrable if and only if (3.5) is also satisfied. 


From Theorem 3.3 and (2.39) follows : 


Corollary 3.3—The distribution D @ {&} is integrable if and only if (3.4) and 
(3.5) are satisfied simultaneously. 
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The existence of maximal elements is proved in closed, bounded and convex, 
but not necessarily compact, subsets of Banach spaces, The theorems proved 
in the paper are general enough to apply to the standard infinite-dimensional 
commodity spaces used in economic analysis, 


1. INTRODUCTION 


Suppose that K is a subset of a Hausdorff topological vector space E. Then each 
binary relation P on K gives rise to a multivalued map 7: K -> 2X as follows: if 
x € K, then 7 (x) = {y € k: (x, y) © P}. Conversely, if T: K > 2* is a multi- 
valued map, then a binary relation P on K is defined as follows: (x, y) € Pif and 
only if y € T(x). A point xo of K is said to be a maximal element of the map 
T: K + 2X, with respect to the binary relation defined above, if T (xo) = ¢. 


Theorems on the existence of maximal elements have important applications in 
mathematical economics. For example, in recent work in general equilibrium theory 
without ordered preferences, the existence of an equilibrium in an abstract economy or 
qualitative game is often proved by constructing a map P, which may be construed 
as a ‘preference map’, on a subset K of a Hausdorff topological vector Space and then 
by showing that there exists a point xo such that P (xo) = 
Aliprantis and Brown!, Gale and Mas-Colell!9, Sonnenschein!8, Borglin and Keiding5, 
Bergstrom?, Schofield!6, Walker2! and Toussaint?°, In this work the assumption is 


made that the preference map is defined on a compact and convex subset of a finite- 
dimensional Cartesian space, Or more generally, 


¢. See, for example, 


a Hausdorff topological vector space. 

Some of these results have been generalized by Yannelis and Prabhakar23 
(Theorem 5.3) who have Proved the existence of max 
convex spaces and by Border4 (Chapter 7) who has prove 
elements in %-compact convex spaces by using the idea of 


imal elements in paracompact 
d the existence of maximal 
an escaping sequence. 

For an extension of some of the 
the reader is referred to the intere 


argument is somewhat removed from 
further here. 


results in the literature to the manifold setting 
Sting paper of Schofield17. Schofield’s line of 
our main theme and will not be considered 
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The object of this paper is to remove altogether the compactness assumptions on 
the domain and codomain of the ‘preference map’. This is achieved by strengthening 
the assumptions on the ‘preference map’. More precisely, we prove that a contracting 
‘preference map’ (definitions follow) which satisfies the other usual conditions, has a 
maximal element in any closed bounded and convex but not necessarily compact subset 
of a Banach space. The theorms that we prove are general enough to apply to the 
standard infinite-dimensional commodity spaces used in economic analysis such as the 
sequence spaces /p (1 < p < ©), the Lebesgue spaces Lp (1 < p < oc) and the space 
M (K) of finite signed Baire measures on a compact metric space K, since these are all 
Banach spaces. 


The economic motivation for the approach used in this paper will be made 
clearer by the following somewhat informal remarks. In R® the existence of maximal 
elements is easily obtained as a consequence of natural economic assumptions (see 
Debru’). The argument goes as follows. Assume that income and all prices are 
strictly positive. Under these conditions the budget set is easily seen to be closed and 
bounded. Now in R*, a closed and bounded set is compact. Hence if preferences 
are continuous (or even upper semicontinuous) in any of the equivalent topologies on 
R®, there is a maximal element in the budget set and the demand correspondence is 
nonempty valued. 


This argument fails in infinite-dimensional spaces because of their greater 
complexity. In a general topological vector space it is not always true that a closed 
and bounded set is compact. In fact, this property holds only for a very special class 
of spaces called semi-Montel spaces (see Wilansky22, p. 90). Since every semi-Montel 
space is semi-reflecxive and a locally convex metric space is semi-reflexive if and only 
if it is reflexive, it follows that every non reflexive Banach space is not a semi-Montel 
space (see Wilansky22, p. 153 and p. 265 and Duffie® for a discussion of semi-reflexive 
and semi-Montel spaces). 


Let us consider the untoward implication of this fact for the space L.. which 
has been widely used in the economics literature. See, for example Bewley®, Brown 
and Lewis® and Toussaint2°, Since Lo. is not reflexive (see Royden!, p. 191) it is 
not a semi-Montel space. Hence the natural argument which shows that a closed and 
bounded set is compact does not extend to it and the existence of maximal elements is 
more difficult to prove (cf. Jones!2). 


The upshot of this discussion is that in infinite-dimensional spaces there is no 
‘natural’ way in which one can prove the existence of maximal elements 


In certain situations the following approach has been used to solve the problem 
Suppose that the commodity space E is the dual of some Banach space G. An example 
is the space L,, which is the dual of Zi. Under these conditions, the Alaoglu- 
Bourbaki theorem (see Wilansky22, p. 130) asserts that a norm bounded set in E is 
relatively compact in the weak * topology. This means that a norm bouded and , 
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closed set in E is weak * compact. Hence, if preferences are continuous, every 
closed and bounded attainable set in E hasa maximal elemnent. This technique of 
obtaining the compactness of closed and bounded sets by the use of a suitable topology 
and then proving the existence of maximal elements has been used by Florenzano’. 


Observe, however, that this approach only works in the case where the Banach 
commodity space has a predual. It does not work for general commodity spaces. 
Conditions for a Banach space to have a predual are given in Holmes!! (pp. 211-14). 


The object of this paper is to suggest and another way in which the existence of 
maximal elements can be obtained. This approach is more general since it does not 
require the Banach space to havea predual. In view of the difficulty of obtaining 
compact sets in infinite-dimensional spaces, it serves a heuristic purpose by suggesting 
that compactness may be needed to prove the existence of maximal elements in budget 
sets or attainable sets if preferences satisfy a contraction condition. 


2. PRELIMINARIES 


We shall use the following notation. If Kis a subset of a Banach space, then 
int K denotes the topological interior of K, co K denotes the convex haul of K and 


co K denotes the closed convex hull of XK, 


Let X be a Banach space and S a bounded subset of XY. Then the Kuratowski 
“measure of non-compactness’ of S, « (S) is defined by 


x (S) = inf {ec > 0: Scan be covered by a finite number of sets with 
diameter no larger than e¢}. 


Observe that if S is compact, x(S) = 0. The proofs of the following two 
theorms can be found in Lloyd!3, [Chapter 6]. 


Theorem (Kuratowski)—Suppose that (An) is a decreasing sequeuce of non-empty 
closed sets in a Banach space such that an = « (An) tends to zero an 7 tends to infinity. 


oo 
Then A = () Anis non-empty and compact. 
n=1 


Theorem (Darbo)—If S is a bounded subset of a Banach Space, then 
«(S)= «(co S). 
Let Y1, Y2 be metric spaces, 


Then a multivalued map Tf} ¥4 "+ 33148 feat to 
be a k-set contraction 


if for all bounded subsets S of ¥;, T (S) is bounded and 
«(7 (S))Sk«(S). The map 7: ¥; > 2%2 


. | is said to be a strict set contraction if it 
is a k-set contraction with k Cae 


Let K be a subset of a topological vector space and T: K > 2 a multivalued 
map. We say that 7 is a preference map if it is generated from a binary relation P on 
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K so that y € T(x) if and onlyif(x,y)E P. A point x in K is said to be a maximal 
element of the preference map Tif T(x) = ¢. 


3. MAXIMAL ELEMENTS 
We now prove the following theorem on the existence of maximal elements. 


Theorem 1—Let E be a Banach space and D a non-empty closed, bounded and 


convex subset of E. Let the preference map P: D + 2° satisfy the following 
conditions: 


(i) for each x € D, P(x) is convex; 
(ii) P is irreflexive, i.e. for each x € Drx & P(x): 


(iii) for each x € D such that P (x) 4 4, there exists y € D such that 
x € int P-1(y); 


(iv) the preference map Pisa strict set contraction, i.e. « (P (S)) < ta (S) 
for all bounded S, where 0 <1 < |, 


Then there exists a maximal element, ie. a point x* € D such that P (x*) = ¢. 


PRooF : Suppose that the theorem is false. Then P (x) ~ ¢ for every x € D. 
As in Martin!4, [Chapter 4] the proof consists of two parts. In the first part, we show 
that condition (iv) implies the existence of a compact convex subset K of D such that 
P isamap from K into 2X. In the second part of the proof, a fixed point theorem is 
applied to the map P. 
Let Ko = Dand forn = 1,2,... Kn = co (P (Kn.1)). We claim that forn = 1, 2, ... 
Kn C Kn-1 and « (Kn) < 1* « (Ko) (1) 


Now Ki C Ko since D is closed and convex. Furthermore, 


a (Ki) = « (co P (Ko)) = « (P (Ko)) <1 « (Ko) 


where the second equality holds because of Darbo’s theorem and the last inequality is 
a consequence of the fact that P is a strict set contraction. Hence, (1) hold form = |. 


Assume that (1) holds forsomen > 1. Then 


Kny1 = co (P (Kn)) C co (P (Kn-1)) = Kn 


and 


a (Kns1) = « (co P (Kn)) = «(P (Kn) < t « (Kn) € t"*1 & (Ko) 


so that (1) holds for alln and the proof of the claim is finished. Since Kn is a 
decreasing sequence of closed sets such that « (Kn) + 0 Kuratowski’s theorems implies 
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oo *. . 
that the set K = () Kn is nonempty and compact. K is also convex as an inter- 
n=0 


section of convex sets. 


Observe that P (Kn) C P (Kn_1) © P (co (Kn-1)) = P (Kn) so that P maps K 
into 2%, and the first part of the proof has been completed. 


It only remains to check that the map P: K — 2X satisfies all the conditions of 
the fixed-point theorm of Tarafdar!®. Now P(x) ~ ¢ for every x € K C D, by 
hypothesis, so that P is nonempty valued on KX. Condition (i) implies that P is 
convex-valued on K. It only remains to prove that for each x € K there exists » € K 
such that x € int P~!(y) in the relative topology of K. So suppose x € K C D. 
Then condition (iii) implies that there exists y € D such that x € int P-1(y). This 
means that y € P(x) C Ksothat y € K. Hence, for each x € K there exists yE K 
such that x € int P-!(y) where the interior is in D. Therefore, there exists an open 
neighbourhood AN in D such that x € N C int P-1(y). Since Nis openin D, N (\ K 
is openin K. Consequently, x € N () K C relative interior of P-1(y) in K. This 
proves that for each x € K there exists y € K such that x € int P~! (y) in the relative 
topology of K. 


Hence, from the fixed-point theorem of Tarafdar!9 we conclude that there exists 
a point x9 € K such that x € P (x9), contradicting condition (ii). The contradiction 
proves the theorem. 


Corollary—Let D be a closed, bounded and convex subset of a Banach space. 
Suppose that the preference map P : D — 2° satisfies the following conditions: 


(i) for each x € D, P (x) is convex; 

(ii) foreach x € D, x & P (x); 
(ili) for each y € D, P~1 (y) is open in D; 
(iv) P is a strict set contraction. 


Then there exists a maximal element. 


PROOF : Itis easily verified that if P-1 (vy) is open in D for y € D, then 
condition (iii) of the the theorem is satisfied. 


A map Satisfying conditions (i), (ii) and (iii) of the corollary is a B-map. The 
following generalization of a B-map is due to Yannelis and Prabhakar23. 


ee Definition—A map T : K — 2K where X is a subset of topological vector space E 
is said to be a BS map if the following conditions hold: 


(i) foreachx € K, x €& co T (x); 
(ii) T~1 (y) is open in K for eachy € K. 
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Clearly every B-map isa BS map. The following theorem is due to Yannelis 
and Prabhakar®8, (Corollary 5.1). 


Theorem 2—Let K bea compact convex subset of a Hausdorff topological vector 
Space and P: K - 2¥ a map that is locally BS-majorized (i.e. locally P is a sub-map 
ofa BS map). Then P has a maximal element. 


The assumption of the compactness of the domain can be weakened in Banach 
spaces if P is a strict set contraction. More precisely, we have the following theorem: 


Theorem 3—Let D bea closed, bounded and convex subset of a Banach spaoe 
E. Let the preference map P: D -> 2? satisfy the following conditions: 
(i) Pisa locally BS-majorised. 
(ii) P is a strict set contraction. 

Then there exists a maximal element. 

Proor : Suppose that the theorem is false. Then P (x) # @ foreach x € D. 
Proceeding as Theorem |, we get a compact convex subset K of D such that 
P:K 2X, Hence, P (x) ~ ¢ for every x © K. Since P is locally BS-majorized, 


for each x, there exists a BS-map Tz ; D — 2° such that z € co 7; (z) for allz € D 
and an open neighbourhood Uz such that z € Uz implies P (z) C Tx (z). Define 


T).(z) = co Tz (z) \ K and U).= Uz; (\ K. Then T..: K + 2* is convex-valued and 


for every z © UL, P(2) CT. Clearly, z & T., (z) for all z € K. Hence, the map 
P: K ~ 2* is locally &S-majorized and Theorem 2 implies that P (x*) = ¢ for some 
x*. The contradiction proves the theorem. 


Remark : Since every locally B-majorized map is locally BS-majorized, it 
follows from Theorem 3 that every locally 8-majorized preference map that is a strict 
set contraction has a maximal element in any closed, bounded and convex subset of a 
Banach space. This enables us to weaken the compactness condition of a Theorem 
2.2 of Toussaint?®, 


We turn now to the consideration of acyclic preference maps. First, we prove 
a generalization of a theorem of Bergstrom? and Walker?!. 


Theorem 4—Let K bea compact subset of a topological space. Suppose the 
preference map P : K — 2X satisfies the following conditions: 


(i)  Pis acyclic; i.e. if 
Xt41 © P (xt) fori = 1,2, .., then xy & P (xn+1); 


(ii) foreach x € K such that P (x) + 4, there exists y € K such x € int P™! (y) 
in K. 


Then there exists a maximal element. 
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Proor : Suppose that the theorem is false. Then P(x) + ¢ for each x € K. 
This implies that for each x € K there exists » € K such that y € P (x). Hence, 
condition (ii) implies that for each x € K there exists y€ K such that.x € int 
P-1(y) = Oy, The relatively open sets O, cover K i.e. K U Oy. Since K is compact, 


n 
there exists a finite subset {y1, yo, ..., ¥n} such that K u Oy s 


It is easily verified that since P is acyclic, the set {1, yo, ..., vn} has a maximal 
element y’ € {)1, 2, -. , Yn}, Le. Jt & P(y’) fori = 1, 2, ...,n. Hence y’ & P-1(y4) 


n 
(a fortiori, »’ & Oy ) for i= 1,2,...,n. This is a contradiction since K = A; Oy 4 
y’ € Kand y’ & Oy for i = 1, 2,...,n. The contradiction proves the theorem. 
{ 


Remark ; If we assume that P~1(y) is open for each y € K, we get the theorem 
of Bergstrom and Walker which is special case of Theorem 4 above. 


For contracting preference maps the assumption of the compactness of the 
domain can be weakened. More precisely we have the following theorem. 


Theorem 5—Let D be a closed bounded and convex subset of a Banach space E. 
Suppose that the preference map P: D > 2” satisfies the following conditions: 


(i) P is acyclic; 

(ii) foreach x € Dsuch P (x) + ¢, there exists 
y € D such that x € int P~1(y) in D; 

(iii) P is a strict set contraction. 


Then there exists a maximal element. 


Proor : Suppose that the theorem is false. Then P (x) ¢ for every x € D. 
Proceeding as in the first part of Theorem |! we get a compact convex set K such that 
P: K — 2X, Arguing as in the second part of Theorem | we can prove that for each 
x € K there exists y € K such that x € int P-1(y) where the interior is in the relative 
topology of K. Hence, condition (ii) of Theorem 4 is satisfied. The restriction of P 
to Kis clearly acyclic. This implies that condition (i) of Theorem 4 is satisfied. 


Consequently, Theorem 4 implies that there exists a point x9 € K such that P (xo) = ¢. 
The contradiction proves the theorem. 
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ON THE ENDL-TYPE GENERALIZATION OF CERTAIN SUMMABILITY 
METHODS 
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It is shown that if A is a regular Euler, Taylor or Meyer-Kénig matrix and 
r> 0, then A and its Endl-type generalization A’ are absolutely equivalent 
for all sequences sy, = o (n'’?), 


Endl4 and Jakimovski® introduced the “Endl-type” generalization of Hausdorff 
and quasi-Hausdorff matrices and, more recently, Kuttner and Parameswaran’ intro- 
duced the generalization of the same type of the Meyer-Kénig — Ramanujan matrix. 
The Hausdorff, quasi-Hausdorff and Meyer-K6énig—Ramanujan matrices, when they 
are conservative, may be considered as “built” around the Euler, Taylorand Meyer- 
Konig methods A = E,, T, and Sy respectively and reflect some of the properties of 
these methods. A similar statement holds good when we consider the Endl-type gene- 
ralizations of all these methods. It is the object of this note to study some properties 
of the Endl-type generalization A’, where A = E,, Tx or Sx. It is shown that A and 
A’ (r > 0) are absolutely equivalent for all sequences sn = o (n1/2) (see below for the 
definitions). 


Definitions—(\) Let O<a <1 andr>oO. The generalized Euler matrix 


E = (of, ) is defined by 
ee 
ay, m= y\ 2 —k) at (1 — ark (k < a) 
0 (k >n) aad) 
n,k = 0,1,...,... when r = O this gives the Euler matrix E,. We write ank for a. 


LetO0<«<1landr>0. The generalized Taylor matrix T7 = (2, jis de- 
fined by 


eee 

Tr - 

ee a n—m an-m (y a mtrty (n ee 0) 
ee 0 (n < m), sas(2) 
*The author thanks the NSERC of Canad 
referee for his remarks on the paper, 





a for partial support of this work, and the 
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When r = 0, this gives the Taylor matrix T,. We write bmn for Das ; 


(3) Let O<a<1 and r>0. The generalized Meyer-Kénig matrix S! 


= ( eg ) is defined by 


m+r+n 


Cnn = en (m +r) = ( ~ 


Ja — ajmtr+1 gn (m,n > 0). 


When r = 0, this gives the Meyer-Kénig matrix S.. We write cmn for ce . {For 


the definitions of Ta, S, given above and further properties of these matrices, see 
Meyer-K6nig!9}. 


(4) Following Cooke?’3 we shall say that two given matrices A and B are absolutely 
equivalent for sequences s of a class X if As — Bs € (co) whenever s € X. {For the 
matrices and sequences we shall discuss in this note, this notion coincides with 
“‘vollaquivalenz” as defined in Zeller and Beekmann15}. Absolute equivalence with 
respect to fixed class XY of sequences is obviously an equivalence relation and is in 
particular transitive. 


Theorem 1—Let0 < « < 1 andr >0.Then E, and E’ are absolutely equi- 


valent for all sequences sn = o (n1/2), 
Proor : Let sn = o(nl/2), Ex, s = t = {tn} and E- s = u = {un}. 


Now 


| tu —un| = | Dank sk — Lai, sk| < ZB | ane — ai, lise |. ...(4) 


n 
(Throughout the proof of Theorem !, the symbol 2 stands for ee .) It is enough to 


prove that the sum in (4) tends to 0 as n — ©9, 


Case 1: r = 1—Since a, Qn,1,k+1, the sum in (4) is 


Z| ank — anzi,ky1\\sk | < B | ank — an,kzillsk | 
= 2 | Qn,k+1 — an..1,k,1\|sk | 
= An + Bn (say). But sk = 0 (k1/2) and hence (i) An + 0 (see Lorentz), 


Theorem 2 and §5.2 or Parameswaran!!, Theorem 3) and (ii) Bn > 0 (Parameswaran??, 
Lemma 5). Hence the theorem is true when r = I. 


Case II: r is a positive integer > |—Now the sum in (4) is 
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r-1 
Z| ank — anyrker|lsk |< UL | anyt,kyt — anytsr,k4t,rllse | 
t=0 


+ 0 ann — co by repeated applications of Case I to {sn} and its translates. 


Case II]: 0 < r < 1—For arbitrary fixed integers n > k > O, let 





ank (x) = ( ae ) okta (| — q)n-k .. (5) 
(n+ x)(n— 1+ x)... (kK + 1+ x) akt7 (1 — an k 

os (n —)! 

Then 
n 
log ank (x) = 7 log (x + v) — log (n — k)! 
+ (k + x) log « + (n — k) log (1 — «) 
and 
‘ n 

Gnk (x) ] 

—_—_ = Sere 4 loge: ... (6) 

ank (x) oS, ee Y ee 


n 
Since log « < 0 ang >> 2 + : is a positive decreasing function of x which tends to 


v~k+] 





+ co asx > — (k + 1) andto0asx > + co, there exists x9 = xo (n, k) such that 
the right-hand side of (6) is positive for —(k + 1)<x < xq and negative for 
x > Xo. Thus 


the function anz (x) is strictly increasing for (*) 
—(k + 1)< x < x9 and strictly decreasing for x > xp. 
ank (t) a(n a t) 


Now, for any t > — k, gap ees 


and this equals | exactly when t = (vn — k)/(1 — a). It follows from (*) that ¢— ] 
< X0 = Xo (n, k) < t; that is, 





on —k —k 

“— I < x0 (n, k) << St med, 
Hence 

Gnk = ank (0) < ank (x) < ank (1) = an ,1,k41 for all x € (0,1) if. < x0 
and hence if 1 < (1 — a)~1 (an — k)— 1. Thus 


0 < ank (x) —~ @€nk < Qn415k41 —Ank 
if 
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; re, ...(8) 
Similarly ank = ank (0) > ank (x) > ank (1) = anyi,t41 for all x € (0, 1) if 0 > xo, 
and hence if 


on — k 
1—e« 





IN 


0. (9) 


Thus if n, k satisfy (8) or (9), then 
| ank (x) — ank | < | an,1,k,1 — ank | . (10) 
for all x € (0, 1). 
The inequality (10) may fail to hold only for those n, k such that neither (8) nor (9) 
holds; thus (10) holds except possibly when 


an — k 


We note that for each n, there are at most two the values of k for which (11) holds. 
We write the sum in (4) as 
x | ank (r) — ank|lse | < 1 + Ze 


where 2; denotes the sum taken over those k which satisfy (10) and 22 denotes the sum 
over those k which satisfy (11). Then 


2X1 S| ank — an,1,k41\|sk | > 0 (asin — cc) as seen in Case I above; and 2 is the 
sum of at most two terms which are 0 (1) since sn = o (n1/2) and max ank (x) = O (n-1/2) 
k 


for each fixed x > 0 (see e. g. Hardy®, Theorem 138 (2) or Lorentz®, p. 313). Hence 
Lo > 0 as n —> ©¢, and the theorem is proved for the case 0 < r < 1 also. 


Case IV—r =[r] + qg, where 1 & [r] = the largest integer less than equal to r 
and 0 < q < 1. Then by (1) 


r (r] 
= | Oo. ~ an, \|sx | 


qd 
mr 2 Qa str) kein] Anstrysk y[ry|lse 


— 0 by Case III applied to the sequence {sx_ry]} (where, as usual, sj is defined to be 0 


when jis a negative integer). Thus Ey and ou are absolutely equivalent for s and 


since E. and ee are absolutely equivalent for s by Case II, we see that the theorem 
is true for the Case IV also. This completes the proof of the theorem. 


Definition 5—We say that a matrix method A is absolutely regular for a sequence 
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s if As — As* €& (co) where s* = sn (n > 0), or equivalently if 


ioe) co 13 
XZ anksk — ZX an,k,1 sk > 0 as n — cc (Cooke?’?). 
k=0 k=0 


Theorem 2—The method E> , where 0 < « < 1 andr > 0, is absolutely regular 
for all sn = o (n1!2) 


Proor : It is enough to prove that 


Pn = Zla,—a’..,, lise | > 0 ifn = 0 (nil), vas BZ) 


This is true when r = O (by Theorems 2 and 9 (ii) of Lorentz®. Since 


Opa a, 2 vor — ank | + | ank — anjky1 | + | anykay 


= Okt [} {se | 
the relation (12) follows from Theorem 1. 
Lemma 1—For 0 < t < l andr = 0, let {E, (t; s)} denote the Ei -transform 
of s = {sn}. Then 
r ?. r *s 
el Uae) be E,4(ss)=t ial (t; 5 — s) 
where 5 = {sn,1}. 


The proof for the case r = 0 has been given in Parameswarn!2; the proof for the 
case r > 0 is similar. 


Theorem 3— Let A = (H’, g) be the generalized Hausdorff matrix defined by 


1 
az = (77) itr (1 — 1)"-* dg (t) 


0 


n 
where g € BY (0, 1] andr > 0, and let s = {sn} = {E ay} satisfy the conditions sn 
0 


= o (nt/2) and an = O(1). Then the set of limit points of u = As will be connected if 
g (1) =¢ (1) =g (1 — 0). 


Remark : Theorem 3 for the case r = 0 was proved by the author in Parame- 
swaran!8; Liu and Rhoades® proved the theorem for bounded {sn} and r > 0; (see also 
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Kuttner and Parameswaran’). Proofs in the more special case sn = O(1) and r = 0 
have been given by various authors; Liu and Rhoades8 or Parameswaran!3 for detailed 


references. 
1 


PROOF OF THEOREM 3: Since un = | ES (t: s) dg (t) we have for n > 1, 


0 


1 1-0 1 
un — Un. = | tE” , (t, a) dg (t) = | 3. | 

: OY eee0 (13) 
o (1) + [g (1) — g (1 — O)J an 


I 


using Lemma | and the facts thnt an = O (1) implies that f E;_, (t, a) is uniformly 
bounded in [0, 1] and tends to 0 as n > co (since sn = o (n1/2) implies that {an} is Ey- 
summable to 0 (Hardy, p. 213) and hence is E. -summable by Theorem 1). It fol- 


lows from (13) that un — un-1 = 0 (1) if g (1) = g (1 — 0) and hence the set of limit 
points of the bounded sequence uv = As is connected, by a theorem due to Barone!. 

Taeorim 4—Let0 <a < land r > 0. Then S, and S’. are absolutely equiva- 
lent for all sequences sn = o (n1/2). 


Proor : The proof is similar to that of Theorem 1 and we omit some of the 


details. Let sn = 0 (n1/2), ¢ = Sy s,u = S. 5,3 = {snyi}, t = SxS and an =n 





haga) eae Ie 
Case I : ris a positive integer—Then um = tmyy by (3) and hence um — tm 
mar ; 
= z (tn — tn-1). Now, as noted by (Meyer-Konig?!9, p. 275), 
N= mt1 
= + 
tn tna =a(l— a > (2F*) ow anys 
v=0 
co 
= ( “—Ja-aen > ("F* ) a Gur — 54) 
itl-—e ahs 


v=0 





I1—a 


-( (tn — tn). But tn — tn). = 0 (1) by the the abso- 
lute regularity of sequences sn = o (1/2) (Parmeswaran', Theorem 3). Since r is a fixed 
positive integer, it follows that um — {tm = 0 (1). 


Case 11:0 <r < 1—It is enough to prove that 
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oo 
x |cn(m +r) — cn (m)ilsn | = 0 (1). 502(14) 
n=0 
We follow Kuttner’s argument used by Sitaraman!4, namely that 
| cn (m + r) — en (m) | < | cn(m) — cn(m + 1) | eee} 
except possibly for those values of » for which 
ma (1 — «)l1 <n < (m + 2) a(l — a)-1. ...(16) 


We write the sum in (14) as £1, + Lo, where 2 is the contribution to the sum by those 

n which satisfy (16) and 2g is the contribution by those n which do not Satisfy (16). 

Now the number g (m) of integers n satisfying (J) isa bounded function of m and 

hence it follows from the well known fact that max cp (x) = O (x-1/2) as x > ©, that 
n 


2) > O0asm > oo, me sf), 
oO co 
Now | Z2| < | cn (m) — cn(m + I)ilsn | = ne cmn — Cm,i,nl|sn | = Pm (say), 
n=-0 n= 
and 
i>?) 
2X cmn — cCm.1, n) sn = 0 (1) whenever sn = o (ni!2) .. (18) 
n=0 


(as we observed in the proof of Case !). But (18) holds if and only Pm = 0(1) Cooke?, 


Theorem 6 or Cooke?, Theorem 5.51). It follows that Xo = o(1). The conclusion (14) 
now follows from (15) and (17). 


Case III: r > 1—The desired result follows from Cases | and II since the mth 
| 
term of S! s is the (m + [r])th term of S ” s, where [r] = the integral part of r. 


Theorem 5—Let0 <a<1andr>0. Then T, and T, are absolutely equi- 
valent for sequences sn = o (n1/2), 


PROOF : We note that the relations (2) and (3) imply that 


” 08 = ~ = . fe = r 
ga = Oforn < mand Om mtn = for all m, n. 
Hence 
co co 
(Or —:Ban)Sq ee (b" —b ) 

mn " m'mtn m,m,n) Smin 

n=0 n=0 
co 
= (co lniace ) 
mn mn) Smin, Fat 4 


n=0 
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We may now prove the theorem either by using arguments similar to those used in the 
proof of Theorem 4, or, as suggested b 


y the referee, we can derive the result from 
Theorem 4 as follows. 


We note that if (amn), (mn) are regular matrices and if {\n} is a given sequence of 
positive numbers then, in order that (amn). (mn) should be absolutely equivalent for 
sequences satisfying sn = o (An), itis necessary and sufficient that 


Rh ehwiee aban FAG =. 0-(1), er) 


n=0 


For we require that the transformation 


oo 
tm = % (umn — Bmn) An (Sn/An) regarded as a transformation from {un} = {sn/ 
n=0 


An} to {tm}, should transform null sequences into null sequences. The “‘standard”’ necess- 
ary and sufficient conditions for this are that (22) should hold and that, for fixed n, 


amn — Bmn — 0 as m — oo, ates) 


But if (umn), (Bmn) are regular, then (23) is necessarily satisfied. Thus, Theorem 4 is 
equivalent to the assertion that 


> |c° = cmp | nl2 = O (1). ...(24) 
mn 
and, by equation (21), Theorem 5 is equivalent to the assertion that 
oo 
SD leg — emn | (m + mi = 0 (1). (25) 
n=0 
Since Theorem 4 has been proved, we know that (24) holds. We have to prove that(25) 


also holds. 


Now choose a constant A with 0 < A < a/(1 — «). Then 


= cmn(m + nil? = O(e-y™) ...(26) 
n< Am 
and 
SS ae (m + n)l/2 = O (e-¥™) eee) 
n<Am 


where y is a positive constant. The first of these results is given by Theorem | = (2) of 
Hardy® the second may be proved in a similar way. But (m + n)1/2 = O (nt?) uni- 
formly inn > Am and hence, by (24), 


706 M. R. PARAMESWARAN 


> | a — cmn | (m+ n)12 = O(1). alae 
n? Am 
The required result (25) now follows from (26), (27) and (28). 
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MATRIX TRANSFORMATIONS IN SOME SEQUENCE SPACES 
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The purpose of this paper is to characterise the matrices in the classes 
(c (p), Cs), (loo (p), Cs) and (1 (p), Cs). 


Let X and Y be any two nonempty subsets of the space of all sequences of com- 
plex numbers and let A = (anx), (n, k = 1, 2, ...) be an infinite matrix of complex 
numbers. We write Ax = (An (x)) if An (x) = © ank xk converges for each n. (Through- 

ke 


out summation without limits runs from | to oo). If x = (xx) € X implies that Ax 
= (An (x)) € Y, we say that A defines a matrix transformation from X into Y and we 
denote it by 4: X¥ — Y. By (X, Y) we mean the class of matrices A such that A:X 
> Y. Ifin X and ¥ there is some notion of limit or sum, then we write (X, Y, P) to 
denote the subset of (¥, Y) which preserves the limit or sum. 


If px > Oand sup pr < o°, we define (see Maddox?) 
k 


c (p) = {x: | xe — 1 | ’e ~ 0 for some /} 


log (p) = {x : sup | xx | % < oo} 
1(p) = (x: % | xe | %e < oo}. 
We define (see Stieglitz and Tietz®) 


n 
Cs ={x: {2 xs} is convergent}. 
i=1 


The purpose of this paper is to characterise the matrices in the classes (c (p), Cs), 
(leg (p), Cs) and (I (p), Cs). 


The following notations are used throughout. For all integers n > | we write 
n 
tn (Ax) = & At (x) = 2 bnk xk 
t=1 


*Present Address : Department of Mathematics, Utkal University, Bhubaneswar 751004. 
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where 


n 
bank = = ak. 
im] 


For all integers n, t > 1 and 1 <r < co, we write 
~ap 


r Ik 
Cin, B, t,r) = Zl bur} B 
ket 


where B is an integer and P, + q, = 1. We put 
C (n, B, 1, ee) = C(n, B) and C (B) = sup C (n, B). 
n 


We have 
Theorem 1—A € (c (p), Cs) if and only if 


-llp 
(i) D= supZ|bne |B * < @ forsome integer B > |! 
n n 


(ii) 3a € C such that 


lim bank = ak(¥ k) 


n7>oo 


(iii) 3 « € C such that 


lm bnk = «., 
n>o 


Proor : Necessity—Let A € (c(p), Cs). Put tn (Ax) = 6n (x). 


Since (c (p), Cs) C (co (p), Cs), {6n (x)} is a sequence of continuous linear functionals 
On co (p) (see Maddox2) such that lim 6n (x) exists. Therefore by uniform boundedness 
principle for 0 < 3 < 1, there exist Ss [0] C co (p) and a constant K such that co» (x) 
S K for each n and x € Ss [0]. Define for each r: 
2 dMlpE sen bank, OS k <r 
f Mars Ovr-ek 


where M = max (1, sup ps). 
k 


Now y? €Ss [0] and 


r —1!/p 
pt lowe | BOY <K 


for each n and r where B == 38 M. Therefore (i) holds. (ii) and (iii) trivially hold. 
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mateo eA? SUpHONS (i)—(iii) hold and x € c(p). Then there exists ] such that 
xh fle <> 0, Hence for 0 < « < 1 there exists ko: ¥ k > ko 
pRiM 
ee eae Sena) < 1 
and therefore for k > ko 


1 
jay 





M |p, MIP 
[xze— I|<B xe — I (55° laa ae dee 
| FOS He 1 res Dal 


By (i) and (ii) we have 


-1/? 
% | buz — 1 | B * = 2D. 


Hence 


= | (bank — ak) (xe — D1 <e. 
k>ko 


Also 


x | (bnk — ak) (xx — 1) | > 0 as n > oc, 
kSko 


Therefore we have 


lim 2 bank xk = ln + & ok (xz — I). 
n>o k k 


This completes the proof. 
Theorem 2—A €E (lq (p), Cs) if and only if 


(i) for allintegers N > | 


iP d ; 
x | bank i N converges uniformly inn 
k 


(ii) lim bnk = ok (¥ &k). 


n7>o 
Proor : Suppose that A € (I. (p), Cs). Clearly (ii) holds. If (i) is false then 


ijt *. . . 
the matrix C = (Cnk) = (ank N *) & (loo, Cs) for some integer V > |, see Stieglitz 
IP, 


1 
and Tietz’. So there exists x € I,. such that Cx & Cs. Nowy = (yk) = (N xk) 
E I, (p), but Ay = Cx & Cs and this contradiction completes the proof. 


For the sufficiency, suppose that the conditions hold. Take an integer N > max 


Vk 
(1, sup | xe | ) We have 


L/pp 
|X (bnk — ak) xKI< aitibanr ee 11 
k 
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By (i) and (ii) we have 


lim 3 bak xk = 2 ak xk. 
n7>o k k 


This completes the proofs. 
Theorem 3—A € (1 (p), Cs) if and only if 
(i) there exists an integer B > 1 such that 


C (B) > co (1 < pr < oo), 

i 
sup | bar | > co (0 < pe € 1) 
nek 


(ii) lim bnz = ox (A k), 
n->co 
This is an immediate consequence of the Corollary to Theorem | of Lascarides and 
Maddox}, if one notices that A € (/(p), Cs) if and only if the matrix B = (bnz) € 
(/ (p), ¢). 
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An exact solution to the unsteady free and forced convection flow of an 
incompressible viscous fluid past an infinite vertical plate is presented and the 
expressions for the velocity, the penetration distance and the skin-friction are 
derived. It is observed that an increase inthe Prandtl number leads to a 
decrease in the penetration distance and the skin-friction when the time f is 
constant, 


1. INTRODUCTION 


Siegell, Schetz and Eichhorn?, Menold and Yang*, Chung and Anderson’, 
Goldstein and Briggs® and Sugawava and Michiyoshi® studied the unsteady free con- 
vection flow under different condition past an infinite vertical plate. Goldstein and 
Eckert? confirmed experimentally some of these theoretical predictions. In all these 
studies, the infinite plate was assumed to be stationary and the fluid was supposed to 
move due to temperature difference only. If the fluid is stationary and the infinite 
plate, surrounded by the stationary fluid, is given an impulsive motion along with its 


temperature raised to Tw such that T|; >T.», where T,., is the temperature of the 


surrounding fluid, how the flow of the fluid takes its shape ? This was studied by 
Soundalgekar® in case of an isothermal plate. The effects of free convection currents 
on the flow and the skin-friction were studied in this paper. However, another physical 
situation which is often experienced in the industrial application is the unsteady free 
and forced convective flow past an infinite vertical isothermal plate of an incompres- 
sible fluid. This has not been studied in the literature. Hence the motivation to under- 
take this study. In section 2, the mathematical analysis is presented and in section 3, 


the conclusions are setout. 
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2. MATHEMATICAL ANALYSIS 


Here we consider the unsteady free and forced convection flow of a viscous 
incompressible fluid past an infinite vertical isothermal plate in the upward direction. 
The x'-axis is taken along the plate in the vertically upward direction and the y-axis 
is taken normal to the plate. Then the physical variables are functions of y’ and t’ 
only. Then under usual Boussinesq’s approximation, the flow is governed by the 
system of equations : 


eu’ ; ; 62 y’ 
Pay = 8 BP (T’ — T) + » Sor (1) 





oT’ é2T 


The initial and boundary conditions are 
t’ & 0, u’ = 0, T’ > T, for all y’ 


fo 0, Ww sO T- aad aty’ = 0. eek) 


---——4- —- YS 


Ue Uo TS = Tak oe oo, 


Here u’ is the velocity of the fluid in the x-direction, P’ the density, g the acceleration 
due to gravity, B the coefficient of volume expansion, 7’ the temperature of the fluid 


near the plate, T, the temperature of the fluid in the free-stream, » the coefficient 


of viscosity, Cp the specific heat at constant pressure and k is the thermal conducti- 
vity. Initially, the plate temperature and the free stream temperature are the same 


everywhere. At t’ > 0, the plate temperature 7). is raised to Te 


On introducing the following non-dimensional quantities 


Y= y UpGriyy ¢ = r’ Us Grlv, u = u’/Uo 


| 
! 
Pr= uCo/k,o = (T'—TLy(r, - 7), | (4) 
| 
ee | 
Gr = vgB Ca Te JUS (the Grashof number) } 
in eqns. (1) — (3), we have 
ou e2u 
or) 7 Ria ...(5) 
Pp 8 tac’ 


or “ay2 ...(6) 
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and the initial and boundary conditions are 


t>0,u=0,6@=Il,aty=0 


t< 0,u = 0,6 = 0, for all yp 
ical Z) 


u=1,@=Oasy—> oo, 


Equations (5) - (6) subject to the conditions ( 7) are solved by the usual Laplace- 
transform technique and the solutions are as follows : 


u = 1 — erfe (4) — 5 ~~ [erfe (n VPr) — erfe (”)] 
+ 2»? (Pr erfe (y Pr) — erfc (x)) 
+ — (exp (—72) — VPr exp (— Pr.n?) ment) 


@ = erfc (y VPr) . (9) 


knowing the velocity field, it is more interesting to study the leading edge effect. The 
penetration distance is derived by integrating u with respect tot and the maximum 
penetration distance Xpmax at any time can be determined by differentiating xp with 
respect to y holding ¢ constant and then by setting the derivative equal to zero. Thus 
the penetration distance is given by 


t 
Xe = Ju (y, t) dt. ...(10) 
0 


This can be expressed in terms of the Laplace transform and inverse transform with 
respect to the variable ¢ as 


L-1 ee Lu, ») } 


Xe 


raf 2 i (y, 5) : eesti) 


Substituting for z (y, s) in (11), and taking the inverse, we have 
2 exp (— 2) 
Xp=t la — erfc (y)) + 7 (ae — 2n erfe (») ) 


oe Sra | Ay2 {(n2 hips + 3 Pr) erfe (nV Pr) 
Ea 


— (72 + 3) erfc (n)} + 3 (erfc (qn V Pr ) — erfe (n)) 


+ 22 ((2n? + 5) exp (— 2) — (2Pr 2 + 5) VPr exp (—Pr 7) | 
os ...(12) 
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To understand the physical meaning of the problem, we have calculated the numerical 
values of xp for different values of Pr and these are plotted in Fig. 1. We observe from 
this Fig. | that the penetration distance xp decreases with increasing the Prandtl 
number of the fluid. But it increases with f, the time. 





Fic. 1 Penetration distance 


We now study the skin-friction. It is given by 





a as | 
py? wae | was) 
and in view of (4), equation (13) reduces to 


du 


2 areg® v'/P US Grillz. 7 
n 








.. (14) 


n=O) 
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Hence substituting for u from (12) in equation (14) and simplifying, we get, 


ib ] ets 1 
— Te ae oy. a ese ou (15) 


From (15), we conclude that the skin-friction decreases with increasing Pr when ¢ is 
constant. 


3. CONCLUSION 


The penetration distance and the skin-friction decrease with increasing the 
Prandtl number. 
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The effect of thermal diffusion on the onset of instability in an unbounded 
vertically stratified, stable, quiescent theromohaline fluid with compensatiug 
horizontal thermal and Salinity gradients saturating a porous medium is 
investigated, Itis found that for this system, instability sets in only through 
stationary mode when the wavelength is large compared to the porosity of the 
medium. The dependance on thermal diffusion of (a) the maximum growth 
rate (b) the ratio of fluxes (c) the onset and region of instability, are shown 
using the Soret parameter, §. For (i) S = 771 — 1, the system remains Stable 
for small perturbations (ii) § < — 1, convection sets in even in the absence 
of the horizontal gradient when both temperature and salinity gradients are 
stable where + js the ratio of mass diffusivity to thermal diffusivity, it is 
also shown that salt gets transported faster than heat when § < 77] — 1 
while transport of heat is more for S> +7! — 1 for some growth rate of in- 
stabiiity. 


1, INTRODUCTION 


The study of thermohaline instability with thermal diffusion in a fluid saturated 
porous medium is of importance in geophysics, ground water hydrology, soil sciences, 
oil extraction, extraction of ores etc., because it is known that the earth’s crust is a 
porous medium saturated by a mixture of different types of fluids like oil, water, gases 
and molten form of ores or ores dissolved in fluids. Thermal gradient present between 
the interior and exterior of the earth’s crust may help convection to set in. Also the 
two transport processes (heat and mass transfer) interfere with each other and produce 
cross phenomena known as thermal diffusion (Soret effect) and diffusion thermo 
(Dufour effect)!’2, Thermal diffusion is the flux of 
dient and diffusion thermo is the flux of heat cause 


"Present Address : Shrimathi Devkunvar 
Madras 600044. 


miass caused by a temperature gra- 
d by a concentration gradient. But, 





Nanalal Bhatt Vaishnav College for Women, 


THERMAL DIFFUSION 717 


in general, in liquid mixtures diffusion thermo is negligible fora heat-solute pair. 


Various authors have analysed Soret driven instability for Newtonian fluid layers?-5 
as well as for fluid saturated porous media§-14, 


All these authors have ignored the horizontal 
nity that may be present. 


of different permeabilities 
lower layers resulting in a 


gradients of temperature and sali- 
The earth’s crust made up of large sedimentary materials 
and heat capacities!5 causes nonuniform heating of the 


aye horizontal temperature gradient. Also, due to the nonuni- 
form distribution of the constituent solutes in the multicomp 
the earth’s crust, 


cal ones. 


onent system saturating 
One encounters horizontal solutal gradients in addition to the verti- 


The aim of the present paper is to study the effects of both thermal diffusion and 
horizontal gradients on thermohaline instability in a fluid-saturated porous medium. 
Holyer!® has analysed double-diffusive (thermohaline) interleaving due to horizontal 
gradients in a Newtonian fluid layer. The present work is based on Holyer’s!® work. 


Numerical computations are made by choosing the values of the Soret parameter 
S at random since it is reported by Hurle and Jakeman®, Schechter et al.4, and Legros 
et al 18, that the sign and magnitude of the Soret coefficient depend on the velocity of 
flow and/or solute concentration of the system. 


2. MATHEMATICAL FORMULATION 


Consider an unbounded region of a porous medium saturated with a Boussinesq 
fluid which is quiescent and which has both temperature and salinity variations. It is 
assumed that the horizontal temperature and salinity gradients compensate each other 
in the basic state and hence there is no change in the density in the horizontal direc- 
tion. The mass flow generated by the temperature gradient is taken into account in 
the formulation of the problem which is analysed for two-dimensional perturbations. 
The coordinate axes are chosen with z-axis vertically upwards and x-axis horizontal in 
the direction of increasing salinity. Only a stable basic state (i.e. where the total 
density gradient is negative) is considered to study the problem. The governing equa- 
tions are 








Ou OW _ eae 
px 1 oz ma 
1 aq See VP Pg 4 

a nit See oo rtd) 
parse, A919 5, Po k 


where A is a constant dependent on the porosity and the geometrical parameters of the 
porous medium. This is the Navier-Stokes equation modified for Boussinesq approxi- 
mation and Darcy law9!?: 


oT 


7a Mq.yT=Ky?T sae 
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ac 
ot 


P = Po[l — ar (T — To) + ae (C — Co)] peal 


49. VC=DC+ DVT (4) 


> 

where u, w are the horizontal and the vertical components of the velocity g: ft, p, », P 

are the time, pressure, kinematic viscosity and density of the double diffusive fluid; ¢ 
> 


and k are the porosity and permeability of the medium; g is the acceleration due to 
gravity; T and C are the temperature and concentration; Po is the density at temper- 
ature Ty and concentration Co; K and D are the thermal and solutal diffusivities; D1 
is the Soret coefficient; xr and ac are the coefficients of thermal and mass expansions 
which are both positive; M = (Po cp)s/(Po cp)m, cv is the heat capacity, suffixes f and 
m stand for fluid and fluid-solid mixtare. But M has been dropped in the subsequent 
analysis because liquid-saturated porous media is considered for this work. 


Small perturbations are imposed on the basic state giver by 


-_> 


qo = 0; To (x, z) = Ta x + Tzz + To; Co (x, 2) = Oext+Qezt+ Co; 
M4) — 
Ps = 0 (by assumption) i. e., ay Tz — ac Cr = O 


where qo, Tp (x, z), Cy (x, z) are the velocity, temperature and concentration in the 
basic state; Tz, 7z are the horizontal and vertical temperature gradients; Cz, Cz are the 
horizontal and vertical salinity gradients. The perturbed state is given by 


Tp (x, z) + T (x, z, t); Cy (x, z) + C (x, z, t); 


Po (z) + p (x, z, t); q (x, 2,1) 


where pp is independent of x from the basic state momentum equation. 





Let the stream function ¥ be defined by u = — a andw = Oy 
az on 4 
The linearised equations in T, C and wv are given by 
F ay = av — 
a 2 hati oa 
( ot Ky ) if Bo Tz ‘eed Tz ...(6) 
7] ‘ 
a7 Pvt)C=p 2 RE Oy SOE: 
ee V wh a res sa ge ..(7) 
ase v aT ac 
ee = 2 a: 
( d ot +- E)Y = @ (ar ax ail? ie ). ...(8) 


On eliminating T and C (6) — (8) lead to 


d ‘ 6 lca 
= 2 ; 
Ga Eve) (tee pve) (2 Waa E)ve 


(equation continued on p. 719) 
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a2 
o x2 


a : 2 
+ gacts( 5° —ky) ox2 





col play 
== gerf, (= =.D.y3) 


+ g ac Cz (K — D) 92 Sez 


+ goc Dj { T, ibe is i 
Ze 1( if: Sa Tz Ths Ey \ve wv. +009) 


The solution to (9) is assumed in the form 
~ = do exp (faz x +iazz +4 1%). 


On substitution for & in (9), the equation in s which decides the stability of the 
System is given by 





(s + Ka2) (cs + Da?) («5+ - ) 


+o ga’ (ar Tz — ac Cz) + ga2 ace Cz [K — D (1 + S)] az az 
— ga®[KacC, — D(1 + S)arT,] a =0 


2 . 


Zz? 


where az and az are the horizontal and vertical wave numbers and a2 = a® +a 
S = ac Di/«r D is the Soret parameter. 


For a?/$ < 1 (i. e., when the wave number is small compared to ¢ or, the wave- 
length is large when compared to ¢. This does not imply that ¢ is large. It is small 
enough to justify Darcy model with ‘a’ < 4), the above equation becomes 


2 a2 
(c + Ka?)(s + Da?) a + aga. (27 Tz — ac C2) 


+ gar%acCz[K — D(1 4+ S)] az az 
— ga?[K ae Cz — D(1 +S) ar Tz] a =0. ...(10) 
For neutral stability through stationary mode (i. e., « = 0) (10) gives 


KD va‘ 
k 


= — gacCz[K —D(1 + S)] az az 
+g [KacC, — D(1+ S)er Tz] a. . mac tid 
From (11) it is seen that ‘a’ will be real when 


we Cr {1 — +t (1 + 5) = eee Cue 2 (1S) tr T's 
aris) 


D 
where +t = K° 
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The inequality (12) is the same for a Newtonian fluid layer (i. e., without porous 
medium) and hence the values of az/az for which convection sets in through stationary 
mode in a Newtonian fluid layer are the same even in the presence of a porous medium. 


For neutral stability through oscillatory mode (c = i, where © real is the fre- 
quency) the equation (10) gives 
wary K Dy a6 


~ —— + ? + garacCz[K — D(1 + S)] az a 


= 4 
— ga2[K ae C, — D(1 + S)or Tz] a tio | = (Kis 


+ h 4 a. (ar T; —- de C2) | = Q. ne (13) 
The imaginary part of (13) equated to zero gives 
4 ms 
@ Eat + Dv +g a. (ar Tz — ae Cz) ] = 0, 


implying » = 0 for a stable basic state (xr Tz — ue Cz > 0). 


Hence, when the wavelength is large compared to ¢ convection does not set in 
through oscilatory mode. Under this assumption, only stationary convection is dis- 
cussed in what follows. Since the Soret perameter S can have any value?*4-18 the dis- 
cussion is based on different ranges of S 


The basic state in which both vertical thermal and salinity gradients are stable 
(ar Tz > 0, we Cz < 0) is considered in the discussion given below. Itis seen from 
(12) that the inequality cannot be satisfied when Cz = Ofor S> — 1. This shows 
that Cz derives the instability and if, in addition, S < t-1 — 1, the inequality (12) 
gives a 0 implying that hot salty solution overlies cold fresh solution which is 

So 

usually the set up for stationary mode. Eut for S <*— 1, convection can set even in 
the absence of Cz provided oc Cz — +(1 + S) ar Tz > 0. A difference in the set up 
with Cz # 0 is noticed for S > t-1 — 1, when the stationary convection sets in with 


> 0 for which cold fresh solution is above hot salty solution which is a diffusive 
x 


regime. This tilting has been reported by McDougall!® for a Newtonian fluid layer. 
It is also seen that (12) is not satisfied for any Cz when S = r-1 — 1, i.e.,a system 
which is stable in the basic state remains stable for small perturbations, implying that 
Cz is compensated by this value of the Soret parameter. A similar analysis can be 
carried out for other directions of the gradients subject to a stable basic state condition. 
The regions of stationary instability in the ae Rt plane are shown in Fig. 1 (a) 
Zz 
— (c) for different values of S and H for given Re 
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Re = -100 


UNSTABLE 


UNS TABLE STABLE 


(b) 


Re — Re =0 Ri- R, =0 


Re = 100 
STABLE 


UNSTABLE 





Biz-ea: 





ao R; plots for « = 0.01 and different values of A. 


(a) —1<§<7+71—1. The region of instability lies between the line considered and 
R: — R, = 0. 
(b) S < — 1. The region of instability is above the line considered. 


(c) S > +7! — 1. Theregion of instability lies between the line considered and R; — R, = 0. 


k 
where the horizontal Rayleigh number H = was ; 
7z k 
the vertical thermal Rayleigh number Ri = el a ae ; 
eS Ge Ge k 


the vertical salinity Rayleigh number Re = — 9 
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It can be seen from these figures that an increase in H increases the region of in- 
stability. Figure | (a) shows that when S increases, the region of instability decreases. 
Figures 1(b)—(c) also reveal that the Soret parameter affects the regions of insta- 
bility. 

3. MaximuM GROwTH RATE OF INSTABILITY 

The growth rate of unstable mode is given by o which has to be real for the in- 

stability to grow via stationary mode. The maximum growth is obtained from 


> 


o 06 
SeaC oa => 0 — ar 
az 0 az 


...(14) 


» 


Proeeeding on the same lines as Holyer!®, the Rayleigh numbers for the maximum 
growth rate are given by 


ae (A+ 1)(A4+ 7) 











Haws aot [EST .. (1Sa) 
(A+ 1) [a2 t(A +1) — 2a” A (A + *)] 
pre ; ...(15b) 
Aa, {lI—+(1 + S)] 
Pe ea : {a2 [A2 (1 — + S) 
Aa [1 —+(1 + S)] (15c) 
+2Atr+22(14 SJ— WAAt+INAt>)} 
where 
pe K a2 


It can be seen that S affects all quantities though Sis independent of S. To 
t 


find the maximum dimensionless growth rate A and the corresponding horizontal and 

vertical wave numbers az, cz, equations (15) are to be slmultaneously solved for given 

and values of the gradients and +, But, following Baines and Gill20 and Holyer16, jt 

is easier to compute H, Ri, Re assigning values for oy and A for given +t. att _ hatin 
. . —~s ax Kt R 

plots a displayed in Fig. 2 (a), (b) and 3 for different values of S for given values of 
Zz 

A and a These figuresare drawn for positive values of Rt subject to ths condition 

Rt — Re > 0 (stable basic state) and for values of Sinthe range —~1<S < 71 

— with +t = 0.01 for a heat-salt pair. Itisseen from Figures 2(a), (b) and 3 that 


wh Ky ° INC P 5 
en S increases R, increases. The effect of S becomes less for higher growth rates, 


For the growth rate \ = 1/2, it is found that ae = | for all <2 and S,i.e ia 
Rt az i eee ie 
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’ 0.10465 


0.0745 





S:1 
OD445 
* 
0.045 
eA 
-0.0455 
0.0120 0.0280 0.0440 0.0640 
t 
28 Re 
Fi. 2. (a) “R- — ~R, plots for R; S @ subject to R;— R, > Ofor + = 0.01, » = 0.02, 
Na a eet tp | 
Re 
Rt 
P H Re ait Ri — Re > for += 0.01, , = 0.09 
Fic, 2 (b) | iw? plots for Ry > 0 subject to Re P T (1 pam! , 


—-1<S<r}— 1. 
- ee plot is the same straight line (not shown in the figure). Figure 3 shows that 
t 


= decreases when | aa | increases and the effect of S is very nearly the same for 
t ag 


. a 
different —. 
az 
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0.05 0.09 | 0.3 037 
R 


H 


a 
Fig. 3. 6 erie 2 plots for Ry > Ojsubject to Ry — Re > O forts = 0.01, —*> 


=-0.1 
az 


— 0.2, — 0.3 and — < § < 171 — J. 
For marginal stability for which A — 0, it is easily seen from (15) that 


H 


Re " 
ae. vo =< = (] + S) with R; -*% ©o 


R 
hence a + 0. When A — 9, it is seen that oe > oe and 7 —> 1. These results 
t 


Zz 
hold good in the absence of porous medium also. 


4. FLUXEs 
For a propagating mode, 


(¥, T, C) = Re [(to, To, Co) exp (i (ae x +iazz +o t)] .. (16) 


where Uo is assumed and real and ¢ is positive for a stationary mode. Rewriting equ- 
ations (3) (for M = 1 and (4) in the form 


oT - = . 
at ft Weis — bs Tz = Kvy?T 


ac 
ar + be Cz — We Cz = DY2C 4 Di v2 T 


and substituting (16), 7 and Co are given by 


io (az Tx = ae Tz) 
gf ans o+ Kar ...(17a) 


(equation continued on p. 725) 
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a se vy 
— D, a® Boe eal : ...(17b) 


The vertical heat flux Fr and vertical salt flux Fc are given by 
az fa pa Ty 


Lg ae o + Ka2 


Ax 5 
az Cz — 2 Vit TF 
Fe = a[ 2Unz ie Ler Cae phe Dy, a (az Tx az za) Ar v2 
o+ Da (sc + Da?)(o + K a2) 
where Fr = < w» T>, Fo = <w C> and <> denotes an average over a wave- 
length. 


The flux ratio is given by 


az 
ae (Re ~ = H(A +4) 


ac Fo 





Pageee Nee 1) = S-o(R, =) 
ag az 


Substituting for R:, Re and H from (15), this becomes 


art ie A+ |) t (18) 
aeFoe08=6CUC AA( I — tS) 4+ abt 


For S < +1 — 1, it is clear that the flux ratio given by (18) is less than unity 
for all values of A and hence salt is transported faster than heat. The flux ratio is 





ay Fr 





=) eB 
Fic. 4. (a) A — E, Pe plots for s = 0,01 for S < 
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greater than unity (i) for all values of A when t-1 -— 1< S Q 17! (ii) for 


A< wnen S > +7! implying that the heat gets transported faster than salt. It 


sea Sw: 

<+S— |] 
4 x 

is further noticed that for S > +1 witha > Saris transports of heat and salt are 


t ° ‘ ° is 
i i irecti = —.——, the flux ratio tends to infinity, i.e., the 
In opposite directions. AtA ome a 


heat flux is very large when compared to the salt flux. Figures (4) a —(c) are A — 


wad = plots in which the above results are shown. Figure 4 (a) shows that the flux 
ae e 


ratio increases when S increases, but decreases when A increases where S < 1-1 — 1. 
For S = t-1 — | there is no convection. It is seen from Fig. 4(b) that the flux ratio 
is greater than unity for all A and that it increases with both S andA when+-1 — 1 








safe by 
ic b> 5:99.5 
} 
1 i 2 3 
F F. 
Fic. 4, (b) 4 — ey plots for + = 0.01 for (c) A= ae plot for > = 0.01 
tTi-{<§< =I. for S>-7). 


<SqQr1l. The discontinuity of the graph at A = a for S > +-1 can be 


noticed from Fig. 4 (c). 


5. CONCLUSION 


This paper describes in detail the results which are due to the presence of the 
horizontal gradient and thermal diffusion in a fluid saturated porous mndium. It is 
shown that when the vertical gradients of both temperature and Salinity are stable (i) 


for S > — 1, there is no convection in the absence of Cz, but for S < — |, condi- 


tional convection sets in. (ii) for—1< § < -7-1~— 1, when Cz derives convection, 


az : : , 
arr. < 0i.e. hot salty solution is above cold fresh solution. (iii) for S > 7-1 — 1, with 
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ie : } az at ee : 
Cz driving stationary convection os. IS positive 1. e. cold fresh solution overlies hot 
zx 


salty solution which is usually a diffusive regime. It is also shown that for S = +-1 
— 1, no convection is possible for all Cz. Also salt gets transported faster than heat 


when S < +71 — | while transport of heat is more for S > +71 for some values of the 
growth rate. 
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Hodograph transformation is employed for steady, plane, viscous, incom- 
pressible constantly inclined twc-phase MFD flows and a partial differential 
equation of second order obtained which is used to find the solution for 
vortex flow. 


1. INTRODUCTION 


Multiphase fluid phenomena are of extreme importance in various fields of 
science and technology such as geophysics, nuclear engineering, chemical engineering 
etc. In recent years, considerable attention has been paid to the study of the multi- 
phase fluid flow system in non-rotating or rotating frames of reference. Multiphase 
fluid systems are concerned with the motion of a liquid or gas containing immiscible 
inert particles. Of all multiphase fluid systems observed in nature, blood flow, flow in 
rocket chamber, dust in gas cooling systems to enhance heat transfer process, move- 
ment of inert particles in atmosphere and sand or other suspended particles in sea 
beaches are the most common examples Naturally, studies of these systems are 
mathematically interesting and physically useful. The presence of particles in a homo- 
geneous fluid makes the dynamical study of flow problems quite complicated. How- 
ever, these problems are usually investigated under various simplifying assumptions. 


Saffman! has formulated the equations of motion of a dusty fluid which is 
represented in terms of large nuniber density N(x,‘ of very small spherical inert 
particles whose volume concentration is small enough to be neglected. It is assumed 
that the density of the dust particles is large when compared with the fluid density so 
that the mass concentration of the particles is an appreciable fraction of unity. In this 
formulation, Saffman also assumed that the individual particles of dust are so small 
that stokes’ law of resistance between the particles and the fluid remains valid. Using 
the model of Saffman, several authors including Michael and Miller, Liu®, Debnath 
and Basu‘ and S. N. Singh ef al.5, have investigated various aspects of hydrodynamics 
and hydromagnetic two-phase fluid flows in non-rotating system. 


Transformation techniques have become some of the powerful methods for solv- 
ing non-linear partial differential equations. Amongst many, the hodograph transfor- 


i 
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mation has gained considerable success on fluid dynamics problems. Ames‘ has given 
an excellent survey of this method together with its application in various other fields. 
Chandna et a/.8 have used the hodograph transformation for steady MFD flows. Also 


Singh et a/.10 have used hodograph transformation in steady rotating MHD flows aad 
obtained some solutions. 


In this paper, hodograph transformation is employed for steady, plane, viscous 
incompressible constantly inclined two-phase MFD flows and a partial differential 
equation of second order obtained which is used to find the solution for vortex flow. 


2. Basic EQUATIONS 


The basic equations of motion governing the steady flow of a dusty, incompres- 
sible, viscous fluid with infinite electrical conductivity in the presence of magnetic 
field are 


diva = 0 a 9 OI 
P [(a@ grad) @] = — gradp+yculHx H+ KN(o—a +ny2a 
(22) 
curl(a7 x H)=0 Bap) 
div (N 2) = 0 ..(2.4) 
m (v, grad) d = K (a — 2) ae) 
divH =0 ...(2.6) 


where @, 2, H, P, P, 1, » are fluid velocity vector, dust velocity vector, magnetic field 
vector, fluid pressure, fluid density, kinematic coefficient of viscosity and magnetic 
permeability respectively; m is the mass of each dust particle, N the number density 
of dust particles and K the Stokes’ resistance coefficient for the particles. 


The situation for which the velocity of fluid and dust particles are everywhere 
parallel, is defined as11 


‘ee ae Edd) 


where « is some scalar satisfying 
a, grad « = 0 .»0(2.8) 
which implies that « is a constant On the fluid streamlines. 


Introducing vorticity function, current density function and Bernoulli function 
defined, respectively, by 
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geek Mas ey Be (2.9) 
px éy 

pare ces (2.10) 
ae Y. oy 

Reap aps (2.11) 


where U2 = us + us , the system of eqns. (2.1) — (2.6) can be replaced by the fol- 
lowing system 








Cur , due 2.12) 
“1, a Lo oe 
ax | oy 
oF OB 
Fy ~ PEua + w OQ Hy — K(w — N)m = — 2B “@a3) 
E ; oB 
v) = — Pw + pQH + K(x—N)u= ay w».(2.14) 
uy He — ug Hy = f (arbitrary constant) PA BY 
” é Cu ra] a 0 a 
WLw(™ gets) tn (me (G-)+ = Z(+))] 
" x ( a 1) a (2.16) 
maT « Cue Cus a( a = ( m )] 
Wy (2 GE + ws oy )+u(m (4) +a 2 wr) 
=K(5- 1 ) us (217) 
oH, CHe 
TS: Cie ke ee.) 


The advantage of this system over the original system is that the order of the partial 
differential equation is decreased. 
We now consider constantly inclined plane flows and let 


non-zero angle between @ and H. The vector and scalar 
diffusion equation (2.15), yield 


«g denote the constant 
products of i and H, using the 


uy Hoe — ue Hy = UH sin oF) =f aye(it, 19) 


ui Hy + uz He = UH cos ap = f cot ag 


Hm a ae, 


where 
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Solving (2.19), we get 


ohh f 
Hy = v2 (C uy — ue), Ho = ve (C ug + uy) -..(2.20) 


Ne C = cot ao is a known constant for a prescribed constantly inclined non-aligned 
ow. 


Using (2.20) in the system of equations (2.9) — (2.18), we have 


“a Gug 
a ay. =— 0 suet) 


so 
1 oe — PE + OL Cue +m) - K(@—N) 


= _— “Ox. ou lame) 
a hal 
pee ee #OL (Cm — ue) + K(a — N) ue = a 


...(2.23) 





maf a 8 uy 0 uy 77 x 
N x (1 ax + ug By )+m (m 5 (4) 
re SCF) Rete 0m 
2 = 5 ) m rem (Pee 


N 
maf a 0 ug a ue ( aC ac. ))] 
abaG ax SRE eisBx i) + ay \N 











ae 
=K( 4 f 1) us tsi) 
2? + wy us) (S44 SB) 4 (Cah — Cut + 2mm 2) 
te uy ue ay 1 oe 
ou )H 2.26 
(3 - eT} <se4 2.209 
dug Om _ 2.27 
ay pe 
3 Cu tm) 2 (Cu rm yo st 2.28 
ral ees eee U2 i ay 


Let the flow variables 1 (x, ¥), ua (x, y) be such that, in the flow region under 
consideration, the Jacobian 


j= ear satisfies 0 << |J| < 
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In such a case, we consider x and y as functions of uy; and ug such that the following 
relations hold true : 











é uy oy eC uy Ox 

ee ptm I cae 

Ox : dug’ Oy 0 uy 
Ot oo ay ee ee (2.29) 

Ox Cu, oy 0 uy 


Employing transformation equation (2.29) in (2 2!) and (2.26), we get 


Ox” eer sslanae) 
Cul ra) uo 


é 7) 
(Cul — Cu, — 2uy ue) (—* -_— ) 


1 6 uy 0 ug 
2 2 ox ie 
+ (uy up + 2Cun us) (S* + i ) ~ 0. .a(2.31) 
The equation of continuity implies the existence of stream function J (x, ¥) so that 
oy ap 
ie = — uo, Oy = U4. : bvel 2a) 


Likewise, equation (2.30) implies the existence of a function L (uw, u2) called the 
Legendre transform of the stream function (x, ¥) such that 


aL aL 
ta geo ® ...(2.33) 


Employing (2.33) in (2.31), we have 











aL 
(uj = us — 2C uy ug) + (2G us — 2C u’ — 4u1 ue) 
eu" 
1 
627 2 2 ez 
Ou Cu ces i Dei ae ) 2 ate 
Ou~ 
2 
...(2,34) 


Now introducing the polar coordinate (U, @) in the hodo 


graph plane i.e. the (ui, v2) 
plane through the relation : 


uy = Ucos 6, ue = Usin 6 


equation (2.34) gets transformed into 


aL 2 #1 1 ez 1 aL 
dU2 UC aUeg. US gtr > ar iat 





2C aL ‘ 
U2 3 


i ae sas ieee 
where @ is the inclination of vector field a. 
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3. Vortex FLow 
A solution of (2.35) is given by 


L = Bz U2 + (A; cos 6 + Bi sin 6) U 


= Bo (u, + uw) )+ Aiur + By us 


where Aj, By and Bo are arbitrary constant and Bs + 0. In this case, 


ol 7) 
uy 


a = — (2B uy + Aj) 


and therefore the velocity field is given by 


eee ete AT _ tie" By 
oe nt ad 5 aie yma 


These relations represent a circulatory flow. 
From (2.20), we get 


ass (¥— Ai? + + Ap 


and 


Ho = —2B2 SIC (« — Bi) — & + Ad) 
a (X — By)? + (y + Ai)? 


The vorticity € and current density Q can be expressed as 


1 
= — = 0, 
E Bo > Q 


- 
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8.1) 


ak ee) 


piel Suoe 


...(3.4) 


3:9) 


From the integrability condition for B with the use of (2.13) and (2.14) and [(3.3)— 


(3.5)], we obtain 


é 
(x — By) 5 (N - a) + (9 + Ai) G (N— 9) + 2(N — «) = 0. 


Solving (3.6), the number density of dust particles N (x, y) is given by 


~ (x — By) (y + Ai) 


where Cj is an arbitrary constant. From equation (2.8) and (3 3), we obtain 


N +a 


a = Co[(x — Bi)? + (y + Ai)" 


where C2 is an arbitrary constant. 


Hence 


Olt Ca f(x — Bi)? + (y + A1)*I. 
weeps Ot AD. 


eel a.0) 


patoct) 


...(3.8) 


...(3.9) 
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Using (3.3) — (3.5) and (3.7) in (2.22) and (2.23) and solving, we get 











RCT ys 
a y , Oaks” pe ———<$<—$—$—$—— C. 
B= F [(x — Bi)? + (vy + A1)?] + Bo In rE es 3 
4B, 
...(3.10) 
where C3 is an arbitrary constant. The pressure P is given by 
se x — B 
P=C3 + (eB)! + OF Asch ee ln ze 
y+ As 
" 

roe PE hw 


In this case the streamlines are given by 
(x — By)? + (y + A1)? = constant 
which are concentric circles. Summing up, we have : 


Theorem 1—If the dust particle is everywhere parallel to the fluid velocity in the 
steady, plane, constantly inclined MFD flow of an incompressible, viscous, two phase 
fluid, then the streamlines are concentric circles and the dust particle number density 
is given by (3.9). Also the velocity, the magnetic field, the vorticity, the current density 
and the pressure are given by (3.3), (3.4), (3.5) and (3.11) respectively. 


4. CONCLUSION 


There are very few exact solutions of two-phase MFD flows. The mathematical 
complexity of the equations governing the flow of an electrically conducting has pro- 
hibited a thorough analysis. To reduce some of the complexity, it becomes necessary 
to make certain assumptions about the inherent properties of the two-phase fluid. 
Furthermore, all the methods of analysis to this date require that we impose some 
restrictions on the angle between velocity field vector and magnetic field vector. In the 
present work, Saffman model for infinitely conducting two-phase fluid flow considering 
constant angle between u and H, called constantly inclined flow, is taken and exact 
solutions of physical importance are obtained applying hodograph transformation. 
Although the scope of the present work is limited, it is believed that by using the 
approach of this paper and the exact solution obtained, work towards boundary value 
problem of practical importance can be pursued. 
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The instability of a heated layer of a viscous fluid confined between two 
horizontal planes and subjected to a variable gravitational field varying 
spatially with height is investigated. For a layer confined between two stress- 
free boundaries, irrespective cf whether gravitational acceleration is increasing 
or decreasing with height, it is shown that (i) the principle of exchange of 
stabilities is valid when the layer is heated from below, and (ii) the layer is 
stable when it is heated from above. In the latter case, the complex growth 
rate of an arbitrary oscillatory mode lies outside of a circle whose radius 
depends on the wavelength of the mode and the Prandtl number ot the fluid 
but not on the Rayleigh number. The underlying characteristic value 
problem is solved approximately for a linearly varying gravity field and stress- 
free boundaries, It is found that in the case of a layer heated from below, 
gravity increasing upward is a destabilizing influence. 


1. INTRODUCTION 


The importance of convection currents in our environment can scarcely be over- 
estimated; we need only look to our immediate environment and note that the circula- 
tion of the Earth’s atmosphere could not be explained without reference to convective 
motions induced by solar heating. However, the idealization of uniform gravity 
assumed in the theoretical investigations, although valid for laboratory purposes, can 
scarcely be justified for large-scale convection phenomena occurring in the atmosphere, 
the ocean or the mantle of the earth. It then becomes imperative to consider gravity 
as a variable quantity varying with distance from a reference point or surface. 


Although Gresho and Sani? have considered a time-varying gravitational field in 
the convection problem of a horizontal fluid layer, the problem of a spatially- varying 
gravitational field in the plane layer case still remains open, Pradhan and Patra4 
considered the problem of onset of thermal instability in a cylindrical shell of self- 
gravitating fluid heated internally and from below i.e. from the inside boundary. In 
case of a shell heated from below, the narrow-gap approximation led to the plane 
layer problem with the gravity field varying linearly with height from the bottom 
surface. This prompted Pradhan and Samal5 to consider the plane layer problem 
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under a spatially-varying gravity field, neglecting viscosity. They found, among other 
things, that if gravity remained downward (upward) throughout the flow domain, 
neutral modes did not exist. A sufficient condition for stability of a layer heated from 
above was that gravity remain directed downward over a sufficiently large part of the 
flow domain. However, a sufficient condition for instability of a layer heated from 
below was that, besides remaining directed downward, the gravity profile must have 
concave curvature throughout the fiow domain. There also existed a circle limiting 
the growth rates of an arbitrary oscillatory mode of the system. These results, which 
are indicative of the behaviour expected in the limit of small viscosity in the fluid, are 
here extended to the case of finite viscosity. 


2. PERTURBATION EQUATIONS 


Consider a layer of an incompressible, viscous fluid statically confined between 
two horizontal boundaries z = 0 and z =d which are maintained at constant 
temperatures To and Tj respectively. The governing equations for flow and tempera- 
ture, under the Boussinesq approximation, are 


Pose +a: vi ) = — W - Pog (z)k (1 —2(T—To)) +uvra 
ab 
= ee, (2) 
and 
Sta yTaeetT 3) 


where P is the density of the fluid, @ = (u, v, w) the fluid velocity, g(z) the gravita- 
tional acceleration in a coordinate system Oxyz having Oz vertical and Ox, Oyina 
horizontal plane, k is a unit vetor in the vertically upward direction; p the pressure in 
the fluid; ». the dynamic viscocity; T the obsolute temperature; « the thermal diffusivity 
and Po the density at temperature 7». 


In the equilibrium state heat transport is by conduction alone and there are no 
velocities in the fluid. The solution of eqns, (1)—(3) is easily obtained in this case as 
the solution of (the basic state is denoted by an over bar). 


es ...(4) 


Vp=8Pok {1 — «(T — To)} ..(5) 


v2 7 =0. ...(6) 
The basic temperature 7 is assumed to depend only on z, so we can integrate (6) 
to obtain 


are | 
T=T 0+ Bz si 
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where 8 = (To — 7})/d is negative, since this is the solution which gives T = To, 
T; atz = 0, d. 


Writing 
p=b+p',T=T+T,P=p+p,a—ia --(8) 


where p’, 7’, P’ and w’ are small perturbations with respect to which the equations are 
linearized and using the Boussinesq approximation, we obtain 


Fi l = 

——_—_—- li : 2 Ti 9 , 

ar po WP tvved+eaeTk (9) 

vy :-# =0 ...(10) 
and 

aT’ 

3 tewarweT .. (11) 


where v is the kinematic viscosity. 


Seeking solutions in terms of normal modes whose dependence on x, y and 1f is 
given by 


exp i (kn x + kyy + ct) wl?) 
and eliminating p’, u’ and v’ from the resulting equations, we have 

yD? )( Do—-)\y  eg ar (13) 
and 

bo TD — 0 (D?.-F 8) 7 = — bw (14) 


where D denotes differentiation with respect to z and k2 = k*. + ke : 
Assuming the boundaries to be perfectly heat conducting and either rigid or free 
we have 
T’ = 0=w for z= 0 andd 
and either Dw = 0 (in case of a rigid boundary) 
or D2 w = 0 (in case of a free boundary) TS) 
Introducing the non-dimensional variables 


ic d2 


24 =2/d, amkd, o, = 
v 


Wy = wdlk, De = dD, 0, = T/pd ...(16) 
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and g (z) = go y (zy) where 80 is the value of g (z) atz = 0, we 


Can write the eans. 
(13) and (14), on suppressing the star subscripts, as : 


(D2 - a2 — oa) (D2 — a’?)w = — Raz (z) @ aah?) 
and 

(D2 — a2 — P,)g=w retin) 
where 

ye —200 Bd" a° (19) 
and 

Py = vjx ...(20) 
are the non-dimensional Rayleigh and Prandtl numbers, respectively. 

The boundary conditions reduce to w = 0 = 6 (z = 0, 1) and either 

Dw = 0 (at a rigid boundary) 
or 

D2 w = 0 (at a free boundary). stel) 


Thus, for a given a, R and Py, the equations (17) and (18) together with the 
boundary conditions (21) constitute an eigenvalue problem for o and the system is 
unstable, neutral or stable according as the real part of o, namely, or is positive, zero 
or negative, respectively. 


3. THE PRINCIPLE OF EXCHANGE OF STABILITIES 


For fixed, stress-free boundaries at z = 0, 1 we have the kinematic boundary 
condition as 


w= D2? w = 0 (z = 0, 1). ed Ay 
Since the thermal boundary condition is 

6=0(z = 0, 1) (23) 
the boundary condition (22) can be expressed in terms of 9 by means of (18) and (23) 
as 

D* 9 = D419=0 (z = 0, 1). pvs( 24) 
So finally 

6 = D29 = D49=0 (z = 0, I). alee) 
Combining eqns. (17) and (18), we have 

(D2 — a2 — P; a) (D2 — a2 — oc) (D2 — a2) 9 = — Ra®y(z)6 _ ...(26) 
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which must be considered together with the boundary conditions (25). 


Multiplying eqn. (26) by 6*, the complex conjugate of 6, and integrating the 
resulting equation over the range of z and making use of the boundary conditions (25), 
we have 


1 1 
Pro { (| DO| 2+ af | | Adz + o (1-H Be) EDS as) | aide 
0 
1 
+ f (| D36@|2+ 3a2| D2 6|2 + 3a4| Do | 2 
0 


1 
+ a8|6|2) dz — Ra® f y(z)|0|2dz =0. Pies) 
0 


Now, for a neutral mode, we must have o = jot with ot real, and then the real and 
imaginary parts of (27) give 


1 al 
Prot § (| Da|2 + a | 6|2)dz — { (| D39|2+4 3a2| D29| 2 
0 0 


1 
+ 3a4 | Da| 2+ a8| | )2dz+ Ra® f y(z)|e|2dz=0 
0 


...(28) 
and 


1 
(1 + Pr) ot | | (D2 — a2)9|2dz =0. ...(29) 
0 


Equation (29) shows that oi = O and then (28) shows that the principle of exchange of 
stabilities is valid provided the neutral state exists, that is, 


1 
J (| D3 6 | 2 + 3a2 | D2 g2 | 2 + 3a4| Do | 2+ a&| 9 | 2) dz 
0 


1 
= Ra? f y(z)|@| 2 dz. ...(30) 
0 


Equation (30) shows that, since y (z) > 0, a necessary condition for the existence of 
neutral states is that 


R> 0. ae ea 


If R is negative, there is no neutral state and a random disturbance will either be 
damped or amplified. 


4. A SUFFICIENT CONDITION FOR STABILITY 


Since c = 0 for R negative in the case of stress-free boundaries, multiplying 
eqn. (27) throughout by 9*, the complex conjugate of o and dividing by | o | 2 we have 
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1 
: 1 
PV OPOVES BLO de + (1 + BT | (De — ayo [2 ae 





1 
o* 
+ rap | (1 D801? + 302] p20) 24 304 | Do] 2 +08 | 0|*) de 
0 


Ra? o* 
—— | y (2) | @|2dz =0. .--(32) 
0 


Separating the real and imaginary parts of eqn. (32), we have 


1 
Pror J (| Do|® + a7 j6|%)dz+(1+P,) | | (D2 —a2)9)2dz 
0 








1 
Eeaeee D3 9 | 24°3@2| D29/2 
|a|2 ( | 6 | + 3a | D 6 | + 3a4 | Do| 2 + a6 | g| 2) dz 
0 
1 
— Ra2 a | +@lel%dz=0 (33) 
0 


and 
1 1 
Prow| (| Del®+a¢\a|%dz—~ —%— [ (| ps9} 
0 0 


+ 3a? | D2 6 | 2 + 3a4| Do | 2 + a6 | @ | 2) az 
1 


a 5 | y¥(z)|6|2dz=0. -+(34) 
’ 0 





+ Ra 


Equation (34) shows that if R < 0, then o, < 0 provided y (z) > 0 over most of the 
layer thus implying stability for the configuration. 


5. A CIRCULAR ExCLUSION THEOREM FOR OSCILLATORY MOopeEs 


Consider an arbitrary oscillatory mode of the system in the case of stress-free 
boundaries. For such modes, we have o¢ # 0. Hence eqn. (34) can be put as 


1 
J (| D3 | 2 + 3a? | D2 6|2 + 3a4| Do|2 + aS|@|2 
0 


u 
eee 1) oF Pr (| D9 | * a8 | 8 | 8) de = 0. 
ED 
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This can be rewritten as 


1 
{ (| D3 612+ 3a2| D20| 2+ 2a4| Dé | 2) dz 
0 
Z 
+ (at = Pr Vo TEC G aes eae 


1 
= Ra2 f y(z)| 6| 2 dz. ..-(36) 
0 
Hence if R < 0, then we must have 
oF + a > allP,. ...(37) 


Hence the complex growth-rate of an arbitrary oscillatory mode (in the case of a 
layer heated from above in a variable gravitational field directed downward) must lie 
outside of a circle whose radius depends on the wavelength of the mode as well as the 
Prandtl number of the fluid (but not on the Rayleigh number of the configuration). 


6. THE PROBLEM OF A LAYER WITH FREE BOUNDARIES 


We solve the characteristic value problem consisting of (17), (18) and (21) for 
free upper and lower boundaries for the case of a linearly varying gravity field 


y¥() =1+Mz>0. ...(38) 


Since the principle of exchange of stabilities is true in this case and the marginal state 
is stationary, the equations to be solved are 


(D2 — a?)2w = — Ra2(1 + Mz)9@ ...(39) 
and 
(D2 — a®)9=w ...(40) 


together with the boundary conditions 

w= D?w = g = 0, (z = 0, 1) : .. (41) 
Eliminating w from eqns. (39) and (40) we get 

(D2 — a?)3 6 = — Ra2(1 + Mz) @. ...(42) 


The kinematic boundary cendition can be expressed in terms of @ by eqn. (24). 
In view of (39) and (40) this gives 


D6 9 = D869 = Dig =... = Dg = 0, (z = 0, 1). ...(43) 


So 6 must satisfy 


= D?¢@= Dig=... = Dg — 0, (z = 0, 1), .+(44) 
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We have to solve eqn. (42) with the boundary condition (44). Rewriting eqn. (42) in 
the manner 


(D? — a?)3 9 = (1 + Mz) v .»(45) 
and 
b= — Rarg ...(46) 
we expand y and @ in the forms 
co co 
a : = — Am sin mz 
U SS Am sin maz and 6 = 33 ay FEES ...(47) 
m=1 m=1 


where @ is the solution of the equation 
6o 
(D2 — a?)3 9 = (1 + Mz) = Am sin m & Zz ~..(48) 
m= 


which satisfies the boundary conditions (44). 


Inserting the expansions (47) in eqn. (46) and multiplying the resulting equation 
by sin nxz and integrating from 0 to 1 givesa set of homogeneous linear algebraic 
equations in the Am. The requirement that the constants Am in the resulting algebraic 
equations are not all zero leads to the secular equation. An approximation to the 
solution can be obtained by setting the determinant formed by the elements in the first 
n-rows and n-columns of the secular determinant equal to zero, successive approxima- 
tions being obtained as n increases. The calculation can be done for various values 
of M, in each case the value of R which corresponds to a particular value of a being 
found by solving an algebraic equation. The minimum value of R with respect to a 
can then be found out for each M. 


A first approximation to the solution of the secular equation is obtained by 
setting the (1, !)-element of the matrix equal to zero and ignoring all the others. This 
corresponds to the choice of sin 7z as a trial function for 6. As Chandrasekhar! has 
shown, very good approximation are obtained for m = | =m and the solution 
converges rapidly for increasing values of mand n. Setting the (1, 1)-element of the 
secular determinant equal to zero, we find that 





b(n? + a2) = =4M+}. Pasi 

This gives 
pa aga CiaSin Bd a ...(50 
R= sy 5 (50) 


We can get a better approximate solution by taking the first two terms in the 
expansion for @ and w. In this case, the secular determinant reduces to the form 
ou Su} (51) 


= 0 
So1 Seo | 
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where 


(72 “+f a2)3 
ve 2Ra2 


8M (822 — a?) 
9n2 (x2 + a?) 


Sie 


So = 


{1 — cosh a + sinh a(l + cosh a)} + 


(4n2 + a2)3 


See 2Ra2 


a2 + 1372 
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iM + 3 
32 (2a4 + 8a2 m2 + 9n4) 
(a2 + 4n2) (a2 + 72)? 


Mn2 (72 + a2) (4n2 _ 3a?) 
a? (4n2 + a2)? 4n2 + a? 


9n2 


4a (1 + cosh a) 
sinh a 


+21M + ft. 


ba RY) 


Equation (51) is quadratic in R. The minimum positive value of R with respect to a 


for each M can be calculated from this 


t,- 


equation. 


DISCUSSION 


Numerical results summarized in Figs. 1—3 were obtained by evaluating (51). 
For comparison purposes, separate curves (dotted ones) have been drawn using (50). 
It is seen that a variable gravitational field with gravity increasing linearly upward has 
a destabiling effect (Fig. |), the neutral stability curves shifting progressively downward 
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Fie. 3 Behaviour of critical wave number a- With increase in gravity parameters. 


as the gravity parameter increases. The monotonic increase in the destabilizing effect 
with increase in the gravity parameter is further brought out in Fig. 2. Figure 3 shows 
that the critical wave-number decreases monotonically with increase in the gravity 
parameter. 


We expect similar behaviour for other boundary conditions. 
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